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I.  Introduction 

This  report  examines  the  operation  of  high-efficiency, 
single-  and  double-dri'ft  GaAs  IMPATT  diodes.  The  mechanisms  responsible 
for  high  efficiency  in  these  devices  have  already  been  identified  and 
include  depletion-layer  width  modulation,1  premature  collection  of 
the  generated  electron  pulse,2  and  movement  of  electrons  at  velocities 
higher  than  the  saturated  velocity  during  a  portion  of  the  RF  cycle.3 
Theoretical  analyses  of  one-sided,  high-efficiency  GaA3  diodes  using 
finite-difference  numerical  algorithms  have  been  published.4  AI30, 
theoretical  analyses  of  double-drift  GaAs  diodes  were  performed,5 
however,  the  restrictive  assumption  was  made  that  the  doping  profile 
on  the  n-side  for  an  optimized  double-drift  structure  is  the  same  as 
for  an  optimized  single-drift  structure.  No  results  have  been 
published  which  present  detailed  simulations  of  the  various  high- 
efficiency  double-drift  structures  and  which  show  the  effects  of 
depletion-layer  width  modulation  and  precollection  on  efficiency 
enhancement  in  these  structures. 

To  simulate  these  diodes,  both  simplified  analyses  and  the 
finite-difference  program  developed  by  Bauhahn  and  Haddad4  are  used 
in  this  report.  By  simulating  many  different  structures  for  different 
RF  drives  and  current  densities,  optimization  criteria  are  developed  for 

P 

GaAs  diodes  operating  near  10  GHz. 

The  simplified  analyses  presented  here  provide  an  inexpensive 
method  of  predicting  the  large-signal  performance  of  the  diodes  which 


1 


agrees  well  with  results  of  the  finite-difference  simulation,  as  far 
as  frequency  of  operation  and  voltage  modulation  levels  are  concerned. 
The  finite-difference  program'  uses  explicit  upwind  drift  and  implicit 
diffusion  current  terms.  It  includes  the  variation  of  the  following 
parameters  with  electric  field:  electron  and  hole  velocities,  electron 
diffusion  coefficient,  and  electron  and  hole  ionization  rates. 

Tunneling  was  not  included  in  these  simulations  and  is  not  considered 
to  be  important  for  the  structures  considered  here.  The  generation 
term  is  explicit  in  the  program.  The  boundary  conditions  are  Dirichlet 
for  the  majority  carrier,  and  for  the  minority  carrier  the  boundary 
condition  is  either  Dirichlet  or  a  mixed  boundary  condition  which 
ensures  that  the  total  current  (drift  plus  diffusion)  is  constant  at 
the  boundary  and  equal  to  a  specified  reverse  saturation  current. 
However,  even  when  the  mixed  boundary  condition  is  selected,  the 
minority  boundary  condition  becomes  Dirichlet  if  the  electric  field 
becomes  negative  at  the  boundary;  in  that  case  the  boundary  concentra¬ 
tion  retains  the  value  that  was  set  when  the  electric  field  was  last 
positive,  until  the  boundary  field  becomes  positive  again. 

As  is  explained  in  Reference  4,  the  program  employs  a  "staggered" 
mesh  with  two  kinds  of  mesh  points,  space-charge  points,  and  field 
points.  This  formulation  has  the  advantage  that  the  finite-difference 
equations  are  conservative  under  electric  field  reversal. 

The  material  parameters  used  for  GaAs  at  300°K  and  500°K" ’ *  are 
listed  in  Table  1.  The  properties  of  the  following  GaAs  IMPATT 
structures  are  presented  in  the  following  order:  uniformly  doped 
single-drift,  single-drift  high-low  and  low-high-low  diodes, 
uniformly  doped  double-drift  structures,  hybrid  double-drift  sturctures, 
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Table  1 

GaAt:  Material  Parameters 
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and  double-Read  double-drift  diodes  (high-low  doping  profile  on  both 
s ides  of  the  junction ). 

2..  Approximate  Analysis 

Approximate  methods  for  analyzing,  high-efficiency  GaAs  I14PATT 

diodes  are  desirable,  especially  since  the  coat  per  solution  of.5  the 

•finite-difference  large-signal  program  can  be  significant  (as  much  as 

$10  for  double-drift, nonpunch-through  structures).  The  simplified 

analysis  is  typically  used  to  determine  approximately  -the  optimum  RE 

drive  for  operation  of  d- .particular  structure -at  a  .particular  frequency; 

then  large-signal  .runs  are  made  and  the  parameters  adjusted  until 

optimum  operation  is  obtained.  The  {simplified5 analysis  could  also  be 

used  to  determine  the  optimum  frequency,  given  a  particular  structure 

and  the  RE  voltage  amplitude. 

In  the  simplified  analysts,,  the  electric  field  profile  Is 

assumed  to  be  piecewise  linear.  The  initial  field  profile  i3  determined 

by  funning  the  finite-difference  dc  program.  Figures  1  and  2  show 

the  field  profiles  included  in  the  analysis.  These  profiles  are 

appropriate  for  the  hybrid  double-drift  case  with  uniform  doping  on 

the  p-side  and. high-low  doping  on  the  n-side.  However,  by  making  the 

slope  of  the  segment  from  x.  =  0  to  x  =  w^  equal  to  the  slope  from 

x  =  w  to  x  =  w  ,  it  is  seen  that  the  analysis  can  be  used  for 
ns  n  • 

structures  with  uniform  doping  on.  both  sides ;  Also,  by  making  w. 

P' 

very  small,  the  analysis  is.  applicable  to  one-sided.,  high-low  or 
uniformly  doped  structures.. 

In  the  approximate  analysis,  it  is  assumed  that  the  field  profile 
moves,  up  and,  down,  but  that  it.  retains  piecewise  iinearity  as  the 
terminal  voltage  is  changed.  Thus  dynamic  space-charge  effects  are 


ignored  (dc  space-charge  effects  are  of  course  included  in  the 
determination  of  the  dc  field  profile).  The  terminal  voltage  is 
given  by 


v(t)  =  Vdc  +  VRF  sin  (ut  +  v)  ,  (1) 

where  t  =  0  corresponds  to  the  half-cycle  point.  At  t  =  0,  an  electron 

appears  at  some  location  x^(t  =  0)  >  0  and  a  hole  appears  at 

x  (t  =  0)  <  0.  (These  positions  are  found  from  previous  runs  of  the 
P 

complete  finite-difference  program,  where  the  positions  of  the 

generated  pulse  maxima  at  half-cycle  were  determined. )  At  t  =  0 

the  peak  field  is  E  ,  or  a  =  1  in  Figs.  1  and  2.  Therefore  the 
*  max 

velocities  of  the  electron  and  hole  at  x  (0)  and  x  (0)  can  be  found 

n  p 

using  the  expressions  in  Table  1  since  the  electric  field  is  known. 

'Hie  particles  are  next  advanced  to  new  positions  according  to: 


x  (At)  =  x  (0)  +  v  (E  )At 
n  n  n  n 

x  (At)  =  x  (0)  +  v  (E  )At 
P  P  P  P 


(2) 


where  the  time  step  At  is  a  small  fraction  of  the  cycle  period  and 
E^jEp  are,  respectively,  the  fields  at  the  electron  and  hole  positions. 

At  t  =  At,  a  new  terminal  voltage  is  evaluated  from  Eq.  1. 

Then  the  value  of  a  must  be  found  such  that  the  integral  of  the  field 
profile  equals  V(t).  The  equations  that  yield  a  for  the  profiles  of 
Figs.  1  and  2  are: 
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At  this  point,  both  the  new  particle  positions  and  new  electric  field 
profiles  are  known  at  the  future  time.  Ikerefore  new  particle 
velocities  can  be  computed  and  the  process  is  repeated.  The  particle 
trajectories  are  computed  until  each  particle  encounters  a  point  where 
the  electric  field  goes  to  aero.  The  time  of  this  event  is  stored  and 
at  the  conclusion  of  the  run  these  two  transit  times  are  printed  out, 
along  with  the  velocities  of  the  depletion  edges  at  the  time  of  collection. 

3.  Single-Drift  GaAs  IMPATTs 

3.1  Uniformly  Doped  Structures .  Extensive  runs  for  uniformly- 

doped,  single-drift  IMPATTs  were  not  carried  out  in  this  study; 

however,  one  case  is  presented  here  which  illustrates  the  operation  of 

these  devices  in  the  precollection  mode. 

Hie  donor  concentration  for  the  device  analyzed  is  1.1*  x  1016  cm-3 

and  the  length  of  the  epitaxial  layer  is  3.2l*  pm.  Figure  3  shows  the 

dc  solution  for  J.  =  1000  A/cra2  and  T  =  500°K;  it  is  seen  that  this 
dc 

device  is  nocpunch  through  under  these  conditions.  The  avalanche 
width,  defined  as  the  point  where  the  hole  current  drops  to  5  percent 
of  its  maximum  value,  is  0,75  pm  for  this  case.  The  depletion  width 
at  dc  is  found  to  be  approximately  2.75  pm. 

A  series  of  large-signal  solutions  was  obtained  for  this  device 
at  differert  frequencies  and  RF  voltage  amplitudes.  The  optimum 
efficiency  was  15.1*  percent,  obtained  for  f  =  10  GHz,  =  971  A/cm2 
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FIG.  3  DC  SOLUTION  FOR  UNIFORMLY  DOPED  IMPATT  FOR 

J.  =  1000  A/cm2  AND  T  =  500°K. 
dc 
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ELECTRIC  FIELD.  1.0E+05  VOLTS  PER  Ch 


and  V__,  =  30  V,  which  is  approximately  one  half  of  the  dc  operating 

nr 

voltage.  The  device  voltage  and  current  waveforms  for  this  solution 
are  shown  in  Fig.  h. 

The  "injected  current"  shown  in  Fig.  Ha  can  be  explained  as 
follows.  If  the  convention  is  adopted  that  positive  electron  and  hole 
currents  flow  to  the  left  while  positive  x  is  to  the  right,  then  the 
particle  continuity  equations  are: 


(0  +  R)  +  - 

q 


3J 

3x 


(9a) 


and 


It  ■  (0  *  ">  - 1  S1 


(9b) 


where  (i,R  are  the  generation  and  recombination  rates  (cm“3-a”1)  and 
q  is  positive  and  equal  to  the  electronic  charge.  Adding  Eqs,  9  and 
integrating  over  the  device  length  w  gives 


<l  f  (of  +  It)  dx  =  2q  f  (n  +  K)  dx  +  •Jp(w) 

0  0 

-  J  (w)  -  J  (0)  +  J  (0)  .  (10) 

n  p  n 

It  is  now  assumed  that  the  device  undergoes  a  periodic  oscillation  of 
period  T  seconds.  Integrating  Eq.  10  from  0  to  T  seconds  and  dividing 
by  T  shows  that  the  term  on  the  left-hand  side  is  zero.  Also,  if 
is  the  average  current  for  the  oscillation,  the  following  is  obtained: 


-11- 


i££sfc&£3tfK!  W8»v v-"* 


UNIFOHML.Y  DOPED 
VHF»30. 

JOO  971.  A/CMh*<2 
FREQUENCY  =  10.0  GHZ 


(b)  PHASE  ANGLE,  DEGREES 

FIG.  4  LARGE-SIGNAL  SOLUTION  FOR  THE  STRUCTURE  OF  FIG. 
WITH  VRF  =  30  V  AND  Jdc  =  971  A/cm2. 
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dc 


-  ifT 
T  Jn 


(Jjx  =  Q,t)  +  Jjx  =  0,t)  ]  dt 

v,t)]  dt  .  (11) 


1  (  T 

-j  [Jn(x  =  w,t)  +  Jp(x 


Substituting  from  Eq,  11  yields 


do 


t  T  r  v 

s  |  (q  (G  +  R)  dx  +  Jp(x  =  w,t)  +  Jn(x  =  0,t) 


dt 

(15) 


If  the  injected  current  is  defined  as 


f  V 

Jinj(t)  -  q  I  (0  +  R)  dx  +  J  (x  =  w,t)  +  Jn(x  =  0,t)  , 

0 


(13) 


it  is  seen  that  one  of  tiie  requirement::  for  having  a  periodic  solution 
is  that  the  time  average  of  the  injected  current  must  be  equal  to  the 
average  terminal  current  J(Jc.  The  reason  for  calling  the  expression 
in  Eq.  13  the  "injected"  current  is  that  it  is  a  measure  of  the  total 
charge  injected  into  the  device;  charge  injected  at  the  boundaries 
appears  in  the  Jp  ,J^  terms  and  charge  appearing  through  generation 
appears  in  the  spatial  integral  of  the  generation-recombination  term. 
For  IMPATTs,  the  generation  term  dominates  in  Eq.  13  and  Jjnj(b)  is 
essentially  the  total  generation  occurring  throughout  the  device  at 
time  t. 

The  current  density  plotted  in  Fig.  4b  is  explained  as  follows. 
At  any  particular  time,  the  total  current  is  constant  at  any  point  in 
the  device  and  is  given  by 


-13- 


&SaJ!S» 


J,  ,Q.  =  J  (*>  +  JJX) 

total  p  n 


+  e  TT 


3E 

at 


(14) 


Integrating  over  the  length  of  the  diode  and  dividing  by  w  yields 


^total 


u  J  lJn(,)  +  VX>1  It 
0 


(15) 


The  second  term  on  the  right-hand  side  is  the  displacement  current  and 
does  not  contribute  to  real  power  flow.  Therefore  the  terminal 
current  plotted  in  Fig.  4b  is  given  by 


“  (  [J  (x)  +  J  (x)]  dx  •  (l6) 

ter  w  j  n  p 

0 

Although  the  displacement  current  term  In  Eq.  15  has  been  removed,  some 
displacement  current  may  appear  in  the  terminal  current  calculated 
according  to  Eq.  16.  w  denoted  the  total  length  of  the  diode  that  is 
simulated;  for  example,  w  is  greater  than  3  pm  for  the  diode  of  Fig.  3. 
However  the  depletion  width  of  the  diode  is  only  2.75  pm  and  furthermore 
the  depletion  width  changes  during  the  RF  cycle.  Hence  the  displacement 
term  of  Eq.  15  does  not  in  general  contain  all  of  the  displacement 
current. 

Figure  5  presents  curves  at  three  points  in  the  RF  cycle  of  the 
electric  field  and  particle  distributions  in  the  device.  At  180 
degrees  in  r'ig.  5a,  the  electron  pulse  has  been  generated  and  is 
drifting  to  the  right.  In  Fig.  5b  at  252  degrees  in  the  cycle,  the 
pulse  maximum  meets  the  depletion  edge  which  is  moving  to  the  left  since 
the  terminal  voltage  in  Fig.  4  is  decreasing.  In  Fig.  5c  at  270 
degrees,  the  electron  pulse  has  been  collected  in  the  undepleted 
region.  At  this  point  of  the  RF  cycle  the  width  of  the  depleted  region 
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is  less  than  2  ym,  considerably  smaller  than  the  dc  value  of  2.75  ym. 

Hence  this  device  is  operating  in  the  precollection  mode. 

Between  270  and  360  degrees  in  the  cycle  the  depletion  edge 
moves  from  a  position  of  1.8  urn  to  2.6  ym.  The  corresponding  time  is 
one  quarter  of  the  period,  or  25  pa.  Therefore  the  average  velocity 
of  carriers  in  the  undepleted  region  during  this  time  period  is 
3.2  x  106  cm/s,  considerably  less  than  the  saturated  velocity  for 
electrons  which  is  5  x  106  cm/s.  Therefore  it  i3  seen  that  predictions 
of  the  performance  of  precollection-mode  devices  based  upon  transit 
times  calculated  using  saturated  velocities  can  be  significantly  in  error. 

A  simple  equation  for  the  movement  of  the  depletion  edge  is 
derived  as  follows.  w,tl^(w)  denotes  the  position  of  the  depiction  edge 
and  the  doping  at  that  position  at  a  particular  time.  Aw  denotes  a 
small  change  in  the  depletion  edge  position,  and  AV  denotes  the  small 
change  in  terminal  voltage  brought  about  by  Aw.  Then 


qnn(w) 

w  x  aE  =  w  x  — -  (Aw)  =  AV  .  (17) 

c 

Therefore  the  rate  of  change  of  w  is 


dw  _ _ e _  dV 

dt  ”  qwNp(w)'  dt 


(18) 


For  example,  at  the  end  of  the  cycle  (cut  =  2u),  using  w  =  2.6  ym, 
llD(w)  =  l.U  x  1016  cm-3  and  dV/dt  =  =  2n  x  1010  x  30,  Eq.  18 

yields  dw/dt  =  3.58  x  106  cm/s.  The  peak  electron  velocity  at 
T  =  500°K,  which  occurs  for  a  field  of  h  kV/cm,  is  1.05  x  107  cm/s. 
According  to  Eq.  18,  to  achieve  the  peak  electron  velocity  in  the 
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undepleted  region  would  require  an  RF  voltage  drive  approximately  2.8 
times  the  value  for  the  solution  of  Fig.  1*,  which  clearly  is  not 
pcrsible  for  this  structure.  However,  Eq.  18  also  indicates  that 
reducing  the  doping  level  in  the  drift  region  would  achieve  higher 
velocities  and,  in  fact,  it  is  shown  later  that  near  peak  electron 
velocities  can  be  obtained  for  high-low  and  low-high-low  doping 
profiles,  which  is  responsible  for  the  high  efficiencies  in  these 
structures. 

Examination  of  the  terminal  current  density  of  Fig.  Ub 
shows  that  the  current  from  l8o  degrees  to  270  degrees  in  the  cycle 
is  small,  because  the  positive  current  induced  by  the  drifting  pulse 
is  canceled  by  negative  current  in  the  undepleted  region  which  is 
necessary  for  the  movement  to  the  left  of  the  depletion  edge.  This 
effect  tends  to  increase  the  efficiency,  since  the  optimum  current 
waveform  would  be  a  spike  at  270  degrees.  However,  the  effect  is 
counteracted  by  the  movement  of  the  depletion  edge  to  the  right  from 
0  to  90  degrees,  which  increases  the  current  during  this  portion  of 
the  cycle.  Therefore  in  uniformly  doped  IMPATTs  the  effect  of 
depletion-edge  modulation  is  largely  capacitive  and  does  not 
significantly  effect  the  efficiency. 

The  approximate  analysis  illustrated  in  Figs.  1  and  2  can  be 
used  to  predict  the  RF  performance  of  this  device.  The  field 
profile  is  chosen  to  match  the  electric  field  resulting  from  the  dc 
solution  shown  in  Fig.  3.  From  the  previous  RF  solution  it  was 
found  that  at  180  degrees  in  the  cycle  the  maximum  of  the  generated 
electron  pulse  is  located  approximately  0.9  x^,  where  x^  is  the 
avalanche  region  width.  The  electron  was  injected  at  this  location  at 
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180  degrees  in  the  approximate  analysis,  the  diode  was  driven  with 
a  10-GHz  signal,  and  various  RF  voltage  amplitudes  were  tried  until 
the  electron  met  the  depletion  edge  with  exactly  a  25-ps  transit  time. 
It  was  found  that  the  proper  transit  time  resulted  when  the  RF 
voltage  amplitude  was  50  percent  of  the  dc  voltage,  almost  exactly 
the  same  as  the  large-signal  result.  The  approximate  analysis 
works  very  well  for  this  case  because,  since  the  doping  is  relatively 
high,  space-charge  effects  brought  ab:  ut  by  khe  generated  charge 
are  minor.  Hence  the  neglect  of  space-charge  effects  in  the 
approximate  analysis  does  not  introduce  significant  error. 

The  approximate  analysis  could  alternately  have  been  applied 
as  follows.  If  it  is  assumed  that  maximum  efficiency  will  result 
from  a  voltage  modulation  of  50  percent,  the  RF  voltage  amplitude 
could  be  set  and  the  frequency  varied  until  the  proper  transit  time 
is  obtained.  The  frequency  resulting  from  this  procedure  would 
again  be  10  GHz. 

3.2  Single-Dri  ft  High- Low  OaAs  IMPATTs .  As  was  seen  in  the 
last  section,  the  doping  level  in  uniformly  doped  GaAs  IMPATTs  is 
too  high  to  obtain  near  peak  velocities  for  electrons  in  the 
undepleted  region.  However  the  drift  region  doping  cannot  be  made 
too  small  because  the  transit  time  would  become  too  large  and  space- 
charge  effects  become  severe.  The  results  of  Bauhahn6  indicate 
that  for  operation  at  10  GHz  the  optimum  drift  region  doping  is 
approximately  5-5  x  10 15  cm-3.  If  Eq.  18  is  evaluated  with  the 
same  parameters  as  before  but  with  Hp(w)  =  5-5  x  1015,  the  velocity 
of  electrons  in  the  undepleted  region  at  360  degrees  in  the  cycle 
is  found  to  be  9*1  x  106  cm/s,  or  just  below  the  peak  velocity. 
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Thus  it  appears  that  =  5.5  x  10 15  cm"3  is  a  good  choice  for 
10-GHz  operation. 

The  first  simulation  of  a  high-low  structure  was  for  a  diode 
previously  analyzed  by  Bauhahn,6  specifically  diode  No.  18  in 
Bauhahn’s  report,  where  an  optimum  large-signal  efficiency  of 
28.9  percent  was  reported.  The  dc  solution  for  the  structure  with 
=  1250  A/cm2  is  3hown  in  Fig.  6a  and  the  doping  profile, 
in  Fig.  6b.  Figure  6c  3hows  an  equivalent  low-high-low  doping  profile 
for  this  diode;  it  is  equivalent  in  the  sense  that  if  a  dc  solution 
is  obtained  for  the  profile  of  Fig.  6c,  the  resulting  avalanche 
width  is  the  same  as  for  Fig.  6b.  Also,  the  electric  field  in  the 
drift  region  is  the  same  for  both  profiles.  The  Qc  given  in  Fig.  6c 
is  the  integral  of  the  doping  in  the  high  region,  in  thi3  case 
0.0b  x  lO-1*  x  6.07  x  1017  cm”2,  Hie  high  region  can  be  fabricated 
with  a  doping  spike,  such  that  the  integral  of  the  doping  within  the 
3pike  equals  Qfl.  For  example,  in  this  report  the  standard  doping 
spike  is  taken  to  have  a  Gaussian  profile  with  0.035  pm  half  width. 

To  realize  the  profile  of  Fig.  6c  using  such  a  spike,  the  height 
would  be  approximately  6.53  x  10 17  cm-3  and  it  would  be  centered  at 
0.22  pm. 

The  choice  of  J,  =  1250  A/cm2  was  made  as  follows.  For 
dc 

highly  efficient  precollection  of  the  electron  pulse  it  is  necessary 
for  the  height  of  the  generated  pulse  to  be  slightly  higher  than  the 
background  doping  at  the  point  of  collection.  When  this  is  true, 
the  electric  field  flattens  out  and  takes  on  a  small  positive  slope  in 
the  vicinity  of  the  pulse.  This  tends  to  bunch  up  the  pulse  as  it 


is  accelerated  into  the  undepleted  region,  producing  a  spike  of 
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FIO.  6  (a)  DC  SOLUTION  WITH  Jd=  =  1250  A/cm2  AND  T  =  500°K,  (b)  DOPING  PROFILE  AND 

(c)  EQUIVALENT  LOW-HIGH-LOW  DOPING  PROFILE  FOR  THE  STRUCTURE  SIMULATED. 
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current  near  270  degrees  in  the  cycle,  exactly  where  it  is  needed 

for  high  efficiency.  Since  the  height  of  the  generated  pulse  is 

directly  proportional  to  J^c,  there  is  a  threshold  dc  current  above 

which  the  height  of  the  pulse  is  greater  than  the  background  doping. 

Figure  7  shows  results  obtained  by  Bauhahn6  indicating  that  this 

threshold  current  is  approximately  1000  A/cm2  for  operation  at 

10  GHz.  If  J,  is  increased  too  much  above  this  threshold  value, 
dc 

the  efficiency  drops  since  the  dc  power  increases  faster  than  the 

generated  RF  power.  Simulations  at  different  values  showed  that 

optimum  efficiency  is  obtained  when  the  diode  is  biased  slightly 

above  threshold,  at  J,  =  1200  A/cm2. 

dc 

Figure  8  shows  a  large-signal  solution  obtained  for  the  doping 
profile  of  Fig.  6b.  The  diode  terminal  voltage  and  injected  current, 
as  defined  by  Eq.  13,  are  plotted  along  with  the  terminal  current 
defined  in  Eq.  16,  The  spike  near  270  degrees  is  brought  about  by 
the  precollection  mechanism  previously  described.  Figure  9a  shows 
the  state  of  the  diode  as  the  electron  pulse  is  accelerated  into  the 
unde pie  ted  region.  The  pulse  relaxes  very  quickly  as  it  is  collected, 
of  the  order  of  the  dielectric  relaxation  time  since  it  enters  a 
region  of  low  field.  This  produces  a  reduction  in  current  as  observed 
in  the  terminal  current  waveform.  As  the  terminal  voltage  recovers, 
the  field  in  the  undepleted  region  increases  such  that  the  electrons 
are  extracted  at  velocities  higher  than  the  saturated  velocity, 
producing  the  second,  broader  peak  in  the  terminal  current.  Figure  9b 
shows  the  diode  state  as  the  electrons  are  being  rapidly  extracted. 

At  this  point  the  electric  field  in  the  undepleted  region  is  approximately 
2  kV/cm  for  which  the  electron  velocity  is  approximately  7.8  x  106  cm/s. 
The  efficiency  resulting  from  this  simulation  was  27.67  percent. 
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FIG.  7  MINIMUM  AND  MAXIMUM  DOPING  AND  MINIMUM  DC  CURRENT 
DENSITY  FOR  THE  PRECOLLECTION  MODE  AT  VARIOUS 


FREQUENCIES. 
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Examination  of  the  terminal  current  of  Fig.  8  shows  that  the 
current  is  rather  large  during  the  first  half  of  the  cycle  which 
reduces  the  efficiency.  This  is  partially  due  to  the  movement  of 
the  depletion  edge  to  the  right  and  is  therefore  capacitive  current. 
However  near  90  degrees,  since  dV/dt  =  0,  the  capacitive  component 
disappears  yet  the  current  is  still  large.  This  can  be  understood 
by  realizing  that  not  all  of  the  electron  pulse  is  extracted  at 
higher  than  the  saturated  velocity;  some  electrons  are  left  behind 
and  are  extracted  at  the  saturated  velocity.  Since  the  transit  time 
at  saturated  velocity  is  long,  all  the  electrons  are  not  collected 
until  well  into  the  first  half  ox'  the  following  cycle.  For  this 
reason,  it  is  next  shown  that  reducing  the  length  of  the  diode 
results  in  higher  efficiency. 

Table  2  presents  results  of  simulations  of  this  structure 
at  various  dc  current  densities  and  RF  voltages.  All  the  runs  are 
for  f  =  10  GHz  and  T  =  500°K. 

The  next  structure  examined  was  the  profile  of  Fig.  10a 
which  is  the  same  as  that  of  Fig.  6b  except  that  the  length  is 
shorter.  The  dc  solution  for  J^c  =  1250  A/cm2  and  T  =  500°K  is  shown 
in  Fig.  10b;  it  is  seen  that  the  device  is  slightly  punched  through 
for  this  shorter  length.  Figure  11  shows  the  terminal  waveforms  for 
the  optimum  large-signal  r®  obtained  for  this  structure.  Comparison 
with  Fig.  8  shows  that,  as  expected,  the  terminal  current  is  smaller 
during  the  first  half  of  the  cycle  when  the  diode  length  is  reduced; 
this  results  in  higher  efficiency.  The  efficiency  calculated  for 
the  run  of  Fig.  11  was  29.35  percent.  Table  3  shows  the  results 
obtained  for  this  structure  at  different  RF  voltages,  all  at  f  =  10  GHz 
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Table  2 


Large-Signal  Results  for  the  Structure  of  Fig.  6b  with 
T  =  500°K  and  f  =  10.0  GHz 


VRF 

hl 

Vdo 

HI 

p 

RF 

(W/cm2) 

n 

(Percent) 

0 

(o/cm2 ) 

B 

(u/cm2 ) 

J, 

dc 

(A/ cm2 ) 

31.15 

1*7.4 

1.40x10" 

16.96 

-28.892 

188.14 

1740 

32.0 

1*3.3 

1.35x10" 

17.93 

-26.408 

194.26 

1740 

33.0 

1*3.6 

7. 76x10 3 

9.97 

-14.261 

198.67 

1780 

29.0 

1*8.0 

l.lOxlO4 

15.61 

-26.257 

193.77 

1470 

30.0 

1*1*.  1 

1.774x10 4 

27.23 

-39.411 

198.84 

1480 

31.0 

1*2.9 

1.277x10" 

20.11 

-26.57 

204.24 

1480 

32.7 

1*3.0 

1.215x10 3 

1.88 

-  2.272 

208.86 

1500 

27.0 

1*7-0 

9.6l3xl')3 

16.52 

-26.373 

200.01 

1240 

28.1 

1*1*.  1* 

1.51x10" 

27.67 

-38.245 

205.98 

1230 

28.5 

1*1*.  1 

1.525x10" 

27.90 

-37.55 

206.48 

1240 

29.0 

1*1*.  0 

1.576x10" 

27.90 

-37.482 

205.47 

1280 

30.0 

1*2.6 

1.138x10" 

21.93 

-25.291 

213.96 

1220 

26.4 

1*1*. 1* 

1.038x10" 

23.34 

-29.775 

212.34 

1000 

27.0 

1*3.8 

1.130x10" 

25.71 

-30.994 

215.34 

1000 

28.0 

1*3.0 

1.096x10" 

25.86 

-27.953 

219.55 

985 

29.5 

1*0.9 

3. 954x10 3 

9.84 

-  9.087 

231.57 

983 

! 
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FIG.  11  URGE-SIG8AL  S0LUTI08  FOR  THE  STRUCTURE  OF 
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Table  3 


Large-Signal  Results  for  the  Structure  of  Fig.  10a 
with  T  =  500°K  and  f  =  10  GHz 


VRF  G 

(V)  (u/cm2 ) 

B 

(u/cm2 ) 

n 

(*) 

Vdc 

(V) 

J. 

dc 

(A/cm2) 

*RF 

(W/em2 ) 

1.0  -143.52 

30.591 

0.12 

1*6.6 

1250 

71.736 

5.0  -119.1*1* 

.166 , 36 

2.55 

1*6.8 

1250 

1.49x10 5 

10.0  -  73.466 

198.1 

6.30 

1*6.9 

121*0 

3.67x10s 

15.0  -  52.937 

205. 11* 

10.21 

1*6.9 

1240 

5.955x10s 

20.0  -  1*1.502 

208.85 

ll*.23 

1*6.9 

1240 

8.3x10s 

25.0  -  35.76 

210.71 

19-30 

1*6.2 

1250 

11.175x10s 

27.5  -  40.553 

211.97 

28.68 

1*1*.  2 

1210 

15.334x10s 

29.8  -  1*0.062 

212.7 

29.22 

1*3.8 

1210 

15.48x10s 

23.0  -  39.561 

213.18 

29.35 

1*3.7 

1210 

15.508x10s 

28.5  -  37.321 

215.03 

28.92 

1*3.3 

1210 

15.157x10s 

29.7  -  5.511* 

236.6 

5.03 

1*0.0 

1210 

2.432x10s 
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and  T  =  500°K.  Figure  12a  shows  a  plot  of  negative  conductance, 
susceptance,  and  efficiency  vs.  RF  voltage;  Fig.  12b  presents  dc 
voltage  and  RF  power  vs.  RF  voltage. 

The  next  diode  studied  is  similar  to  diode  No.  12  in 
Bauhahn's  report,6  for  which  an  optimum  large-signal  efficiency  of 
38.3  percent  was  reported.  Figure  13a  shows  the  low-high- low  doping 
profile  for  this  structure  and  Fig.  13b  presents  an  equivalent 
high-low  piofile.  If  a  Gaussian  doping  spike  of  0.035  um  half 
width  were  used  to  realize  the  highly  doped  region  of  Fig.  13b,  it 
would  be  centered  at  0.39  um  with  a  height  of  3.6?  x  1017  cm-3.  The 
high  doping  of  Fig.  13a  is  les3  than  the  value  given  in  Bauhahn's 
report;  this  reduction  was  necessary  to  obtain  the  3ame  electric 
field  profile  reported  there.  The  reason  for  this  is  that  different 
boundary  conditions  were  used  for  these  simulations  than  for  the 
simulations  reported  by  Bauhahri.  Bauhohn  used  Dirichlet  boundary 
conditions  for  both  carrier  concentrations;  the  simulations  in  this 
report  U3e  Dirichlet  boundary  conditions  for  the  majority  carrier, 
however,  for  the  minority  carrier  a  current  boundary  condition  is 
imposed  such  that  the  current  at  the  boundary  is  equal  to  a  specified 
reverse  saturation  current,  at;  explained  in  the  introduction  of  this 
report. 

Figure  13c  shows  the  dc  solution  for  J^c  =  1000  A/cm2  and 
T  *  500°K  obtained  with  the  doping  profile  of  Fig.  13a.  It  is  seen 
that  the  diode  is  significantly  punched  through  at  dc.  An  important 
parameter  for  these  structures  is  EtQ,  defined  as  the  electric  field 
value  at  dc  at  the  beginning  of  the  drift  region  on  the  n-side. 

For  Fig.  13c,  E  =  2.67  x  10s  V/cm.  For  the  previous  structure  of 
Fig.  10b,  E  was  1.8  x  10s  V/cm, 

uO 
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OF  (a)  WITH  T  =  500°K  AND  J  =  1000  A/cm2 


Figure  lU  shows  the  device  waveforms  for  the  most  efficient 
(n  =  35.12  percent)  solution  obtained  with  f  =  8.5  GHz.  The  efficiency 
at  10  GHz  is  not  as  good  because  the  device  is  too  long  for  optimum 
operation  at  10  GHz.  The  peak  in  the  terminal  current  near  l8o 
degrees  is  induced  by  the  generated  electron  pulse  drifting  away 
from  the  avalanche  region;  at  180  degrees  the  device  is  punched 
through.  The  dip  in  terminal  current  near  230  degrees  is  caused 
by  the  inward  movement  of  the  depletion  edge,  when  the  device  becomes 
nonpunch  through.  Near  270  degrees  the  inward  movement  of  the  depletion 
edge  ceases  and  the  pulse  is  accelerated  into  the  undepleted  region, 
both  effects  causing  a  peak  in  the  terminal  current.  As  before,  when 
the  electric  field  recovers  in  the  undepleted  region,  the  electrons 
are  extracted  at  faster  than  the  saturated  velocity.  At  the  end  of 
the  cycle  the  diode  is  punched  through  again;  however,  a  significant 
current  still  flows  because  not  all  the  electrons  have  been  extracted. 

If  this  structure  is  operated  at  10  GHz  the  current  at  the  end  of  the 
cycle  increases,  which  reduces  the  efficiency. 

Table  shows  the  results  of  the  simulations  of  the  doping 
profile  in  Fig.  13a  at  f  =  8.5  GHz  for  different  RF  voltages.  To 
obtain  optimum  efficiency  at  10  GHz  the  diode  was  shortened,  resulting 
in  the  profile  of  Fig.  15a.  The  dc  solution  with  .1^  =  1200  A/cm2 
is  shown  in  Fig.  15b.  Figure  l6  shown  terminal  waveforms  for  the 
most  efficient  (n  =  3**.l8  percent)  solution  at  f  =  10  GHz  and  Table  5 
shows  simulation  results  for  different  RF  voltages.  Figure  17a  shows 
a  plot  of  negative  conductance,  susceptance  and  efficiency  vs.  RF 
voltage  for  this  structure;  Fig.  17b  shows  dc  voltage  and  RF  power 
vs.  RF  voltage.  Comparison  with  Fig.  12  for  the  slightly  punched 
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Table  *4 


Large-Signal  Results  for  the  Structure  of  Fig.  13a 
with  f  =  8.5  GHz  and  T  =  500°K 


v 

RF  G 

(V)  (u/cm2) 

B 

(li  /cm2 ) 

n 

ill 

Vdc 

(V) 

J, 

do 

(A/cm2 ) 

p 

RF 

(W/cm2 ) 

5.0  -73.726 

23.018 

1.00 

77 .2 

1190 

921.51 

10.0  -79. 57*4 

105.18 

*4.27 

77.1 

1210 

3.9785x10’ 

15-0  -63.358 

13*4.14*1 

7.6*4 

76.9 

1210 

7.1276x10’ 

20.0  -51.762 

1*16.16 

10.91 

76.7 

12*40 

10.352x10’ 

25.0  -Itl. 912 

15*4.39 

1*4.39 

76.1 

1200 

13.097x10’ 

30.0  -35.95 

157.62 

17.72 

75- *4 

1210 

16.178x10’ 

35.0  -31.18*4 

159.63 

21.15 

7*4.3 

1220 

19.10x10’ 

I4O.O  -27.51 

160.7 

2*4. 7*i 

72.7 

1220 

22.008x10’ 

*45.0  -23.551* 

162.28 

28.5*4 

70.  <4 

1190 

23.8*48x10’ 

50.0  -21.9**3 

161 . 87 

33.98 

67.3 

1200 

27**428x10’ 

55.0  -18.7*4 

1*47.21 

35.12 

70.3 

1150 

28.3*43x10’ 

SOLUTION 


CURRENT  DENSITY  TERMINRL  VOLTAGE 

K  J  LOflMPS/SQ-CM  VOLTS 


PRECQLLECTION  MODE 
VRF=51.9 

JDC=  1160.  A/CMh*2 
FREQUENCY  =  10.0  GHZ 


PHASE  ANGLE,  DEGREES 


a 


PHASE  ANGLE,  DEGREES 

FIG.  1 6  LARGE-SIGNAL  SOLUTION  FOR  THE  STRUCTURE  OF  FIG.  15a 
RESULTING  IN  n  =  3*4.18  PERCENT. 
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Table  5 


Large-Signal  Results  for  the  Structure  of  Fig.  15a 
with  f  =  10  GHz,  T  =  500°K  and  E  =  2.68  x  105  V/cm 

wO 


v 

RF  G 

B 

n 

Vdc 

Jdc 

p 

RF 

(V)  (u/cm*) 

(u/cmJ) 

(V 

(V) 

(A/cmJ ) 

(W/cm* ) 

1.0  -125.84 

79.127 

.07 

77.3 

1200 

62.898 

5.0  -10U.05 

15*1.96 

1.1*1 

77.2 

1200 

1.3x10 3 

10.0  -  70.1*5 

183.3 

3.76 

77.2 

1210 

3. 522x10 3 

15.0  -  51.58 

193.3 

6.27 

77.0 

1200 

5. 8026x10 3 

20.0  -  1*0.985 

196.63 

8.78 

76.8 

1210 

8. 197x10 3 

25.0  -  3*1.05 

198.2 

11.35 

76.5 

1230 

10.61*xl03 

30.0  -  28.863 

199.2 

13.98 

75.9 

1220 

13.0xl03 

35.0  -  2** . 386 

200.0 

16.63 

7*1.9 

1200 

ll*.  937x10 3 

1*0.0  -  21.73 

200.0 

19.5>* 

73.6 

1210 

17.381*x103 

1*5 .0  -  18.997 

200.0 

22.18 

71.1 

1220 

19.23l*xl03 

50.0  -  16.686 

200.92 

25.71 

68.0 

1190 

20. 857x10 3 

51.0  -  16.729 

201.0? 

27.21* 

67.1 

1190 

21. 756x10 3 

51.8  -  17.01* 

201.23 

29.05 

66. 5 

1180 

22. 862x10  3 

51.9  -  19.1*2 

199.92 

3U .  18 

66.1 

ll60 

26. 155x10 3 

52.0  -  19.121 

200.19 

33.79 

66.0 

ll60 

25. 851x10 3 

52.8  21*. 69 

220.63 

-51.07 

56.2 

1200  • 

-31*.  1*16x10  3 
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through  case  with  smaller  E  shows  that  as  the  diode  becomes  more 
punched  through,  the  efficiency  increases,  the  dc  voltage  rectifi¬ 
cation  becomes  larger,  and  the  RF  voltage  range  for  which  high 
efficiency  is  obtained  is  reduced. 

Previously  published  results7  indicated  that  GaAs  IMPATTs 

should  not  be  punched  through  at  any  point  of  the  cycle.  This 

argument  is  based  on  the  fact  that,  for  a  given  electric  field 

modulation  in  the  avalanche  region,  the  maximum  RF  voltage  swing 

is  obtained  for  an  under-punched-through  structure;  also  the  back-bias 

effect  is  reduced  if  nonpunch-through  diodes  are  used.  The 

simulations  in  this  report  consistently  show  that  back  bias  is 

indeed  more  significant  for  punched-through  structure;  and  that  the 

field  modulation  in  the  avalanche  region  is  greater.  However, 

allowing  larger  field  modulations  results  in  greater  calculated 

efficiencies.  That  this  is  true  can  be  explained  as  follows.  To 

achieve  the  precollection  .node,  according  to  the  simplified  analysis 

(ignoring  space-charge  effects)  the  electric  field  profile  of  Fig.  15b 

must  drop  at  270  degrees  in  the  cycle  so  that  the  field  i3  zero  at  the 

location  of  the  pulse  maximum  at  that  time.  If  E.  is  raised  but 
v  to 

the  field  in  the  avalanche  region  remains  the  same,  the  amount  by 
whicn  the  field  must  drop  from  the  dc  value  at  270  degrees  is 
increased.  This  produces  an  increase  in  RF  voltage  swing.  However, 
raising  E  also  increases  the  dc  voltage  which  increases  the  dc 

wO 

power  (J^c  remains  the  same  for  precollection).  However,  the  dc 
voltage  increase  is  due  to  the  increase  of  the  electric  field  in  the 
drift  region  only,  since  the  field  in  the  avalanche  region  stays 
the  same;  in  contrast,  the  RF  voltage  increase  involves  the  entire 
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electric  field  profile,  since  the  entire  profile  must  drop  to  the 

level  required  for  precollection.  Hence  as  the  diode  is  made  more 

punched  through,  the  RF  voltage  increases  faster  than  the  dc 

voltage,  which  increases  the  efficiency.  This  effect  becomes  more 

pronounced  for  double-drift  GaAs  IMPATTs,  since  then  the  field  on  the 

p-side  remains  unchanged  as  E  is  increased.  (Of  course  the  back- 

bias  effect  must  also  be  taken  into  account  which  reduces  both  VD„ 

nr 

and  V^c  £is  predicted  from  a  simplified  analysis.) 

Another  reason  efficiency  increases  as  E  increases  is  related 

LO 

to  the  length  of  the  device.  Obviously  any  pushed- through  structure 
is  shorter  than  the  corresponding  nonpurched- through  structure.  Ac 
previously  mentioned,  some  electrons  are  left  behind  in  the  pre¬ 
collection  process  and  must  be  extracted  at  saturated  velocity.  If 
the  device  is  too  long,  these  electrons  produce  a  significant  current 
at  the  beginning  of  the  next  cycle,  thus  reducing  the  efficiency. 

The  punched-through  structure  is  shorter,  reducing  the  current  at 
the  beginning  of  the  cycle  and  increasing  the  efficiency. 

The  increased  back-bias  effect  in  punched-through  IMPATTs 
can  result  in  thermal  instabilities,  'ibis  problem  in  addressed  later 
in  the  report. 

Because  the  range  of  RF  voltage  for  which  high  efficiencies 
are  obtained  is  quite  narrow,  as  shown  in  Fig.  17a,  an  attempt  was 
next  made  to  modify  the  structure  of  Fig.  15a  in  order  to  stabilize 
the  mode  over  a  wider  RF  voltage  range.  After  collection  of  the 
pulse  near  270  degrees  in  the  cycle,  it  is  desirable  that  the  depletion 
edge  move  at  near  peak  velocity  for  as  long  a  time  as  possible  during 
the  last  quarter  cycle.  If  wq  denotes  the  position  of  the  depletion 
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edge  upon  collection  of  the  pulse  and  tc»  the  time  of  collection,  then 
it  is  desirable  to  have 

w(t)  =  +  v(t  -  tc)  ,  t  =  tQ  ,  (19) 

where  v  denotes  a  velocity  greater  than  the  saturated  velocity  and  less 
than  the  peak  electron  velocity  and  w  denotes  the  depletion-edge 
position.  Since  the  diode  waveform  is  known  from  previous  simulation 
of  the  structure  with  a  uniformly  doped  drift  region  and  it  is 
a33umed  that  the  voltage  waveform  will  not  significantly  change  when 
the  structure  is  modified,  Eq.  18  may  be  solved  for  the  H^x) 
necessary  to  yield  the  desired  depletion-edge  movement.  Since  dV/dt  =  0 
at  2J0  degrees,  the  doping  must  be  made  smaller  at  the  point  of  collection 
in  order  to  increase  the  depletion- edge  velocity  soon  after  270 
degrees  in  the  cycle.  Hear  360  degrees  the  dV/dt  term  approaches  its 
maximum  value  30  ND(x)  must  reach  a  larger  value  near  the  right-hand 
contact. 

Figure  l8a  shows  a  doping  profile  resulting  from  the  solution 
of  Eq.  18;  it  is  similar  to  the  profile  of  Fig.  15a  except  that  the 
drift  region  doping  has  been  tailored  for  more  efficient  extraction 
of  electrons.  Figure  lGb  shows  the  dc  solution  for  J^c  =  1200  A/cm2. 
Figures  19a  arid  b  show  large-signal  results  vs.  P.F  voltage.  It  is 
seen  that  tailoring  the  drift  region  doping  has  resulted  in  high- 
efficiency  precollection  for  a  wider  range  of  RF  voltages.  Also, 
the  maximum  efficiency  is  slightly  higher.  Hence  there  is  an  advantage 
to  designing  high-efficiency  GaAs  IMPATTs  with  a  positive  doping 
gradient  in  the  drift  region.  Figure  20  shows  terminal  waveforms 
for  this  structure  which  yield  the  best  efficiency,  n  =  35  percent. 


FIG.  18  (a)  DOPING  PROFILE  TAILORED  FOR  MORE  EFFICIENT  EXTRACTION 

OF  ELECTRONS  AND  M  DC  SOLUTION  OF  THE  PROFILE  WITH 

j  =  1200  A/C»l.  T  =  500°K  AND  E^  =  2-67  *  105  »/«• 
dc 
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LARGE-SIGNAL  SOLUTION  FOR  THE  STRUCTURE  OF  FIG.  l8a 
WITH  T  =  500°K  RESULTING  IN  35.0-PERCENT  EFFICIENCY. 
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The  next  structure  studied  was  similar  to  Fig.  18a  except 
that  the  avalanche  region  width  was  reduced  by  one-half.  Figure  21a 
shows  the  resulting  profile  where  the  height  of  the  highly  doped 
region  was  also  adjusted  so  that  the  field  in  the  drift  region  is  the 
same  as  for  the  structure  of  Fig.  l8a.  (it  is  necessary  to  increase 
this  doping  because,  since  the  avalanche  region  width  is  smaller, 
the  maximum  electric  field  in  the  device  is  higher.)  Figure  21b 
shows  the  dc  solution  for  this  profile  with  J  =  980  A/cm2.  If 
the  highly  doped  region  were  realised  with  a  Gaussian  doping  spike 
of  0.035  wm  half  width,  the  height  of  the  spike  would  be  5*5  x  1017  cm’ 
and  it  would  be  located  at  0.2  pm.  Table  6  shows  the  results  of 
large-signal  simulation::  of  this  structure;  it  is  seen  that  halving 
the  avalanche  width  has  increased  the  maximum  efficiency  from  35 
to  38.3  percent.  Figure  22  shows  plots  of  the  large-signal  results 
vs.  RF  voltage  for  f  =  9.5  GHz  and  T  =  500°K. 

The  next  structure  examined  was  a  modification  of  the  structure 
of  Fig.  15a  to  show  the  effect  of  increasing  the  doping  in  the  drift 
region.  Figure  23a  shows  the  profile  with  Hj  in  the  drift  region 
increased  to  7.8  x  10 15  cm-3;  the  device  length  is  also  increased  so 
that  it  is  nonpunch  through  during  the  entire  RF  cycle.  Figure  23b 
shows  the  dc  solution  for  this  structure.  The  maximum  large-signal 
efficiency  obtained  was  22  percent  since  the  generated  electron 
pulse  is  no  longer  higher  than  the  background  doping  at  the  point  of 
collection;  as  previously  explained,  the  pulse  must  be  higher  than 
the  background  doping  for  efficient  precollection  to  occur.  This 
result  confirms  that  the  optimum  doping  level  in  the  drift  region  for 
precollection  at  10  GHz  is  approximately  5  x  1015  cm-3. 
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Table  6 


Large-Signal  Results  for  the  Diode  of  Fig,  21a 


at 

T  =  500°K  and  f 

=  9.5 

GHz 

v 

RF 

G 

B 

n 

Vdo 

Jdc 

p 

RF 

(V) 

(u/cm  *) 

(y/cmJ ) 

(*) 

(v) 

(A/cmO 

(W/cm*) 

5.0 

-92.214 

163.97 

1.64 

72.1 

977 

1.1526x10s 

10 

-61.892 

183.29 

4.37 

72.0 

984 

3.0945x10’ 

20 

-35.895 

191.7 

10.05 

71.7 

996 

7.1788x10’ 

30 

-24.11*7 

193.25 

15.84 

70.8 

968 

10.866x10’ 

40 

-19.098 

193.15 

22.72 

68.3 

985 

15.279x10’ 

1+6 

-16.231 

195.06 

27.04 

65.7 

967 

17.172x10’ 

1+8 

-15.566 

195.89 

28.92 

64.4 

962 

17.932x10’ 

50 

-15.397 

196.65 

31.83 

63.0 

961 

19.246x10’ 

50.5 

-15.685 

196.76 

33.48 

62.5 

956 

20.0x10’ 

51.0 

-17.44 

196-24 

38.3 

61.9 

957 

22.68x10’ 

52 

-11.775 

200.52 

27.17 

59.8 

979 

15.919x10’ 
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LARGE-SIGNAL  RESULTS  FOR  THE  DIODE  OF  FIG.  21a 


(a)  DOPING  PROFILE  MODIFIED  FROM  FIG.  15a  BY  INCREASING  OVERALL  LENGTH  AND  DR. 
REGION  DOPING  AND  Cb)  DC  SOLUTION  WITH  T  =  500°K  AND  Jdc  *  1250  A/cm2. 


h_.  Double-Dri  ft  GaAs  IMPATTs 

^.1  Uniformly  Doped  Structures.  The  approximate  analysis 
developed  in  Section  2  can  be  used  to  design  uniformly  doped, 
double-drift  GaAs  IMPATTs.  From  the  simulations  of  single-drift, 
uniformly  doped  diodes  it  can  be  inferred  that  precollection  of  .the. 
electron  pulse  will  occur  on  the  n-side  near  270  degrees  in  the 
cycle.  Therefore  a  reasonable  design  procedure  is  to  ensure  that  the 
injected  electron  in  the  approximate  analysis  has  about  a  one-quarter 
period  transit  time  before  it  meets  the  depletion  edge.  Also,  it  has 
generally  been  found  in  these  simulations  that  optimum  efficiency 
occurs  when  V^,  is  approximately  one-half  V^;  this  value  of  voltage 
modulation  may  also  be  incorporated  in  the  approximate  analysis. 

Figure  2h  chows  an  '-lec.trie  field  profile  which  resulted  from 

the  approximate  analysis  for  operation  at  10  GHz.  From  dc  solutions  of 

similar  structures,  it  was  found  that  the  peak  field  is  approximately 

U . 0U  x  105  V/cm  at  500°K.  Also,  it  is  a  typical  design  procedure 

to  make  the  electric  field  on  the  left-hand  side  approximately  0.25  E 

max 

this  ensures  that  the  device  is  punched  through  on  the  p-side  for 
most  or  all  of  the  RF  cycle.  The  lengths  in  Fig.  2h  are  obtained 
from  the  approximate  analysis.  If  a  10-GHz  sinusoidal  voltage  of 
magnitude  V^/2  is  impressed  across  the  device,  it  is  found  that  a 
dc  depletion  width  of  approximately  2.1  pm  on  the  n-side  yields  an 
electron  transit  time  of  25  ps.  As  explained  in  Section  2,  the  electron 
transit  time  is  the  time  between  the  electron  injection  at  the  half¬ 
cycle  point  and  the  time  when  the  injected  electron  meets  the 
depletion  edge  on  the  right-hand  side,  where  the  electric  field  goes 
to  zero  in  the  approximate  model.  A  transit  time  of  25  ps  is  desired 
for  a  10-GHz  signiil  because  the  electron  pulse  should  be  collected 


near  270  degrees  in  the  cycle  for  optimum  operation.  Of  course,  at 
the  time  of  collection  the  depletion  width  on  the  n-side  is  signifi¬ 
cantly  smaller  than  it  is  for  the  dc  solution.  (The  initial  electron 
position  at  180  degrees  was  determined  from  previous  large-signal 
simulations  of  similar  structures;  it  was  0.2  pm  to  the  right  of  the 
junction.  The  hole  starting  position  was  1.4  pm  to  the  left  of  the 
junction,  which  was  also  determined  from  large-signal  simulations.) 
The  length  of  the  p-side  was  chosen  so  that  the  hole  transit  time  is 
45  ps,  slightly  less  than  half  the  period.  The  efficiency  is 
improved  if  the  hole  transit  time  is  somewhat  less  than  T/2,  since 
then  the  pulse  is  more  thoroughly  collected  by  the  end  of  the  cycle. 
(Due  to  back  diffusion,  some  of  the  holes  arrive  after  the  pulse 
maximum  at  the  left-hand  boundary.)  An  additional  micron  of 
undepleted  material  is  added  on  the  n-uide  30  that  the  diode  is  never 
punched  through  on  the  n-side. 

A  doping  profile  that  yields  a  dc  electric  field  similar  to 

that  of  Fig.  24  for  J,  =  1000  A/em?  is  shown  in  Fig.  25a;  Fig.  25b 

dc 

shows  the  dc  solution  for  this  structure  at  T  =  500°K.  Tab.e  7 
shows  the  results  of  large-signal  simulations  of  this  structure;  the 
optimum  frequency  was  found  to  be  9  0Hz  rather  than  10  GHz,  although 
better  than  15-percent  efficiency  was  obtained  at  10  GHz.  It  is 
believed  that  these  results  represent  realistic  efficiencies 
obtainable  from  uniformly  doped  GaAs  double-drift  diodes  at  X-band. 

Figure  26  shows  the  terminal  waveforms  for  the  optimum 
solution  obtained  with  this  structure  at  f  =  9  GHz.  Near  180  degrees 
both  pulses  drift  toward  their  respective  contacts,  producing  a 
positive  induced  current  as  shown  in  Fig.  27a;  however,  this  positive 
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Table  7 


Large-Signal  Results  for  the  Double-Drift  Diode  of  Fig.  25a 

with  T  =  500°K 

V  V  P  J 

f  RF  dc  RF  n  0  B  dc 


(Site) 

ill 

ill 

(kW/cm2 ) 

(Percent) 

(u/cm2 ) 

(u  /cm2) 

(A/cm 

9.0 

70.0 

116 

18.8 

16.07 

-7.671* 

98.ll* 

1010 

10.0 

70.0 

117 

17.733 

15.^7 

-7.238 

112.88 

981* 

10.0 

65.0 

119 

18.118 

15.1(3 

-8.577 

111.  ll» 

986 

10.0 

6o.o 

122 

17.568 

ll».  76 

-9.76 

109.7 

979 
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26  LARGE-SIGNAL  SOLUTION  FOR  THE  DOUBLE-DRIFT  DIODE 
OF  FIG.  25a  AT  T  =  500°K  YIELDING  16.07-PERCENT 
EFFICIENCY. 
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FIG.  27  DIODE  STATE  AT  TWO  TIMES  DURHG  THE  BF  CTCIi  OF  FIG. 


current  is  offset  by  a  negative  current  brought  about  by  the  inward 
movement  of  the  depletion  edge  on  the  n-side.  Hence  the  terminal 
current  near  l8o  degrees  is  not  very  large.  Near  270  degrees  the 
electron  pulse  is  collected  in  the  undepleted  region  while  the  hole 
pulse  continues  drifting  to  the  left  as  shown  in  Fig.  27b.  The 
large  terminal  current  between  270  and  360  degrees  in  Fig.  2 6  is 
brought  about  by  both  the  outward  movement  of  the  depletion  edge 
on  the  n-3ide  and  the  movement  of  holes  onthep-side.  The  current 
due  to  the  depletion-edge  movement  has  a  large  capacitive  component 
but  it  is  not  entirely  capacitive  since  the  previously  collected 
electron  pul3e  is  also  being  extracted.  As  with  the  one-sided 
uniformly  doped  GaAa  structure,  the  doping  on  the  n-side  of  double¬ 
drift,  uniformly  doped  diodes  is  too  large  to  achieve  greater  than 
saturated  velocities. 

4.2  Hybrid  Double-Drift  Structures.  This  section  considers 
double-drift  OoAs  IMPATTs  which  are  uniformly  doped  on  the  p-slde 
with  a  high-low  (or  low-high-low)  profile  on  the  n-side.  It  is  not 
possible  to  simply  append  a  uniformly  doped  p-side  onto  a  doping 
profile  that  has  been  optimized  for  a  cae-sided  GaAs  IMPATT.  This 
is  because  adding  the  p-side  significantly  widens  the  avalanche 
region,  thereby  lowering  the  peak  electric  field  at  the  junction; 
hence,  the  electric  field  profile  on  the  n-side  is  quite  different 
from  the  one-sided  case  if  the  same  profile  is  used  in  a  hybrid 
structure. 

The  approximate  analysis  developed  in  Section  2  can  be  used  to 
design  hybrid  structures.  First,  however,  the  dc  program  must  be  run 
at  the  expected  temperature  and  dc  current  density  of  operation  in 
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order  to  approximately  determine  the  maximum  electric  field.  Also, 

an  estimate  must  be  made  of  the  positions  of  the  electron  and  hole 

pulse  maxima  at  l8o  degrees  in  the  RF  cycle  either  through  experience 

running  the  large-signal  program  for  similar  structures  in  the  frequency 

range  of  interest  or  from  an  analytical  expression  for  the  injection 

angle.  Further,  the  value  of  E,  ,  or  the  electric  field  at  the 

to 

beginning  of  the  drift  region,  must  be  assumed ;  initially  it  will 

be  taken  as  approximately  1.8  x  10G  V/cm.  The  procedure  is  then 

to  optimize  the  structure  such  that  the  proper  transit  times  are 

obtained  for  both  carriers.  On  the  n-side  the  electron  pulse  is 

precc llected  due  to  the  inward  movement  of  the  depletion  edge,  much 

like  the  mode  of  operation  in  one-sided  precollection  diodes. 

Therefore  the  electron  transit  time  on  the  n-side  should  be 

approximately  0.25  T.  (A3  previously,  the  electron  transit  time  in 

the  approximate  analysis  is  defined  a n  the  time  interval  between 

injection  at  half-cycle  and  the  instant  the  injected  electron  meets 

the  depletion  edge  on  the  right-hand  side.)  Onthep-side,  the 

hole  pulse  i3  collected  in  the  normal  IMPATT  mode;  optimum  efficiency 

is  obtained  when  the  hole  transit  time  is  approximately  0.U  T. 

Figure  28a  shows  a  hybrid  doping  profile  designed  to  be 

nonpunch  through  during  the  entire  RF  cycle;  Fig.  28b  shows  the  dc 

solution  for  J,  =  1500  A/cra2  and  T  =  500°K.  It  will  be  recalled 
dc 

that  for  the  one-sided  case  the  threshold  dc  current  density  for 
precollection  at  10  GHz  was  approximately  1000  A/cm2.  However, 
for  the  hybrid  structure  the  avalanche  region  width  is  considerably 
wider  and  hence  the  generated  electron  pulse  is  wider.  Therefore 
a  higher  J^c  is  necessary  so  that  the  height  of  the  electron  pulse 
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28  (a)  HYBRID  DOPING  PROFILE  DESIGNED  FOR  10-GHz  OPERATION  AND  (b)  DC  SOLUTION 

AT  T  =  500°K  AND  J,  =  1500  A/cm2. 


is  greater  than  the  background,  doping  at  the  point  of  precollection, 
which  raises  the  field  in  front  of  the  pulse  and  results  in  more 
efficient  collection  of  the  pulse.  The  broader  generated  pulses  and 
higher  required  J^,  as  well  as  the  presence  of  the  p-side  where 
precollection  does  not  occur,  reduce  the  efficiencies  obtained  from 
hybrid  structures  below  those  obtained  with  one-sided  precollection 
mode  IMPATTs . 

If  the  maximum  dc  electric  field  profile  from  running  the  dc 

program  is  approximately  known  (and  the  hole  and  electron  positions 

at  180  degrees  from  large-3ignal  run3  are  known),  the  approximate 

analysis  can  be  used  to  determine  the  dimensions  of  Fig.  28a  for 

10-GHz  operation  as  follows.  A  10-GHz  RF  voltage  waveform  of 

amplitude  V^/2  is  impressed  across  the  device  and  the  lengths  are 

adjusted  until  the  proper  transit  times  are  obtained  for  the 

injected  hole  and  electron,  i.e.,  approximately  40  ps  and  25  ps, 

respectively.  The  electric  field  on  the  left-hand  side  is  set  to 

0.25  E  during  the  adjustments.  It  was  <*ound  that  the  dimensions 
max 

of  Fig.  28a  yielded  a  hole  transit  time  of  36  pB  and  an  electron 
transit  time  of  25  ps. 

Figure  29  3hows  terminal  waveforms  for  the  most  efficient  RF 

solution  obtained  at  10  GHz  with  n  =  18.72  percent.  The  V_„/V 

«r  dc 

ratio  for  this  solution  was  0.6,  whereas  0.5  optimum  voltage  modulation 
was  assumed  in  the  approximate  analysis.  It  is  seen  that  there  is 
a  significant  terminal  current  flowing  during  the  first  half-cycle. 

To  see  why  this  is  the  case,  Fig.  30  shows  the  electric  field  and 
particle  concentrations  at  two  points  in  the  cycle,  306  degrees  and 
0  degrees  (or  360  degrees).  In  Fig.  30a  the  electric  field  reaches  a 
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FIG.  29  LARGE-SIGNAL  SOLUTION  FOR  THE  DIODE  OF  FIG.  28a 
AT  T  =  500°K  WITH  n  =  18.72  PERCENT,  VDT,  =  62.5 

nr 

J.  =  1U70  A/,  n2 ,  V  =  10U  V,  G  =  -  lU . 73  t5/cm2 
dc  ac 

B  =  98.12  u/cm2  AND  f  =  10  GHz. 
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small  constant  value  in  the  undepleted  region  at  approximately 
b. 75  um;  in  Fig.  30b  the  depletion  edge  is  at  7  pm.  The  time 
interval  corresponding  to  this  depletion  edge  movement  is  (360  -  306) 
x  T/360,  or  15.0  ps.  Thei’efore  the  average  velocity  of  the  depletion 
edge  during  this  time  interval  is  1.50  x  10  7  cm/s,  which  is  50  percent 
greater  than  the  peak  electron  velocity.  At  306  degrees  the  electric 
field  in  the  undepleted  region  is  2.3*«  kV/cm  which  corresponds  to  an 
electron  velocity  of  8.67  x  106  cm/s.  Hence  the  "knee"  of  the  electric 
field  is  moving  considerably  faster  than  the  electrons.  For  this 
reason  many  electrons  are  left  behind  in  a  high  field  region  where 
they  are  extracted  at  saturated  velocity.  This  effect  is  brought 
about  by  the  increased  space  charge  in  the  device.  Since  the 
threshold  dc  current  density  for  efficient  precollection  is  now 
approximately  1500  A/cm?  rather  than  1000  A/cm2  for  the  previous 
one-sided  case,  the  width  of  the  depleted  region  on  the  n-side  is 
increased;  therefore,  the  depletion  edge  must  move  faster  to  cover 
a  longer  distance  in  the  came  time  period.  One  way  to  reduce  the 
current  during  the  first  half-cycle  is  to  reduce  E  for  the  dc 

wO 

solution  of  Fig.  28b  which  is  accomplished  by  increasing  the  doping 

in  the  "high"  region  on  the  n-side  in  Fig.  28a,  as  in  the  doping 

profile  of  Fig.  31a.  Figure  31b  shows  the  dc  solution  with 

J.  =  1500  A/cm2  and  T  =  500°K  for  this  case.  It  is  seen  that  the 
dc 

depletion  edge  is  now  at  6.3  pm  instead  of  'i.O  pm  for  the  previous 
case.  Figure  32  shows  the  most  efficient  solution  obtained  for 
this  structure  at  3.0  GHz.  Although  the  terminal  current  during  the 
first  half-cycle  is  somewhat  reduced,  the  efficiency  is  reduced 
because  the  voltage  modulation  is  now  only  0.1*77 .  The  fact  that 
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VRp/V^c  decreases  with  E  is  intuitively  expected,  as  previously 

discussed  for  the  one-sided  case;  this  is  because  reducing  E  reduces 

xo 

the  RF  swing  of  the  entire  electric  field  profile,  but  the  dc  field 
profile  is  only  affected  in  the  drift  region  on  the  n-side.  There¬ 
fore  the  RF  voltage  decreases  faster  than  the  dc  voltage,  as  is  seen 
by  comparing  Figs.  29  and  32. 

To  increase  the  efficiency,  the  diode  of  Fig.  28a  can  simply 

be  made  shorter,  and  hence  punched  through,  so  that  the  transit  time 

is  shorter  for  any  carriers  that  are  extracted  at  saturated  velocity. 

In  Fig.  33a,  the  length  has  been  reduced  to  6.0  pm;  Fig.  33b  shows 

the  dc  solution  for  this  structure  with  J,  =  1500  A/cm2  and 

dc 

=  117.6  V.  The  optimum  RF  solution  for  this  structure  at 

10  GHz  resulted  in  n  =  22.1(9  percent,  V  =  6l.O  V,  =  101  V, 

J^c  =  lU^O  A/cm2,  0  =  -  17.8,  and  B  =  106.5  ti/cm2. 

To  further  increase  the  effic'  y,  EtQ  may  be  increased, 

which  is  accomplished  by  reducing  the  magnitude  of  the  highly  doped 

region  on  the  n-side  in  Pig.  33a.  However,  it  is  expected  that 

increasing  E  will  result  in  a  greater  back-bias  effect,  as  was 
xo 

the  case  for  the  one-sided  structures.  Figure  3*(a  shows  a  doping 
profile  designs  so  that  E  =  2.23  x  105  V/cm;  an  equivalent  Gaussian 

xo 

doping  spike  for  the  highly  doped  region  with  half-width  0.035  pm 
would  have  a  maximum  value  of  U . 05  x  1017  and  would  be  located  about 
x  =  3.1  Pm.  Figure  3*«b  shows  the  dc  solution  with  J^c  -  1500  A/cm2 
and  V^c  =  130  V.  Table  8  presents  the  results  of  RF  simulations  of 
this  structure  at  T  =  500°K,  and  Fig.  35  shows  the  most  efficient 
solution  obtained  at  10  GHz.  As  expected,  the  V^p/V^c  ratio  is  larger 
than  before,  equal  to  0.72.  Also  the  current  in  the  first  hail-cycle 
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Table  8 


Large-Signal  Results  for  the  Hybrid  Structure  of  Fig.  34a 
with  T  =  500°K 


v 

RF 

G 

B 

n 

Vdc 

J, 

dc 

p 

RF 

f 

(V) 

(y/cm2 ) 

(y/cm2 ) 

(1) 

ill 

(A/cm2 ) 

(W/cm* ) 

(GHz) 

10 

-55.621 

62.701 

1.1*3 

129 

1500 

2.7806x10s 

10 

20 

-44.376 

83.96 

1* .  6l 

128 

1500 

8.871*5x10’ 

10 

30 

-34.897 

9I4. 1*81* 

8.24 

127 

1510 

15.703x10s 

10 

40 

-28.0147 

100.13 

11.98 

121* 

1510 

22.1*36x10’ 

10 

50 

-23.1*07 

103.35 

15.86 

122 

1520 

29.257x10’ 

10 

60 

-19.624 

105.86 

19.76 

118 

1510 

35.321x10’ 

10 

70 

-16.155 

107.66 

23.52 

112 

1500 

39.577x10’ 

10 

76 

-lli.921 

108.93 

27.61 

106 

11*50 

1*3.089x10’ 

10 

77 

— lU . 587 

109.23 

27.67 

107 

11*60 

1*3,21*1x10’ 

10 

78 

-14.301 

109.1*9 

27.16 

107 

1500 

1*3.501x10’ 

10 

79 

-13.821* 

109.96 

26.90 

107 

1500 

1*3.135x10’ 

10 

80 

-13.61(5 

110.0 

26.59 

106 

151*0 

1*3.661x10’ 

10 

82 

-12.118 

111.65 

25.51 

106 

1510 

1*0.739x10’ 

10 

81t 

-10.31*7 

113.38 

23.36 

105 

11+90 

36.502x10’ 

10 

86 

-  9.0035 

111*.  36 

20.57 

lOi* 

1550 

33.293x10’ 

10 

88 

-  5.7761* 

117.36 

ll*.  38 

103 

1510 

22.365x10’ 

10 

89 

- 

■263.1* 

77 

-15.11* 

100.98 

27.75 

107 

1510 

1*1* .  88x10’ 

9.5 

77 

-15.091* 

115.93 

26.98 

107 

15l*0 

1*4.742x10’ 

10.5 

77 

-15.07!* 

93.233 

27.18 

108 

1530 

44.684x10’ 

9 
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FIG.  35  LARGE-SIGNA*  ..OLUTION  FOR  THE  HYBRID  STRUCTURE 

OF  FIG.  34a  AT  T  =  500°K  YIELDING  n  =  27.67  PERCENT. 
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is  significantly  reduced  due  to  the  shorter  length  of  the  device. 

Figure  36  shows  the  results  for  f  *  10  GHz  and  T  =  500°K  plotted  as 
a  function  of  Vnr,.  The  negative  conductance  curve  of  Fig.  36a 

nf 

decreases  monotonically ,  whereas  the  negative  conductance  curve  of 

Fig.  12a  has  a  local  maximum  near  the  maximum  efficiency  point. 

Therefore  these  results  indicate  that  it  would  be  easier  to  operate 

the  double-drift  diode  as  an  oscillator  at  the  maximum  efficiency 

point  than  the  single-drift  diode.  The  monotonically  decreasing 

character  of  the  conductance  for  the  double-drift  diode  is  due  to  the 

presence  of  the  p-layer.  In  Fig.  36b  it  is  seen  that  the  dc  operating 

voltage  drops  by  almost  20  percent  from  the  dc  value.  This  in  worth 

noting  since  the  device  heats  up  considerably  for  3mall  VRR,  where 

the  efficiency  is  much  less  and  the  dc  power  i3  increased. 

In  the  doping  profile  of  Fig.  37a,  the  doping  In  the  highly 

doped  region  is  further  reduced  30  that  E  is  increased  to  2.68  x  105  V/cm. 

x>o 

The  equivalent  Gaussian  doping  spike  for  this  cane  would  have  a 
maximum  value  of  3.16?  x  1017  cm-3.  Figure  37b  presents  the  dc 
solution  for  -  1500  A/cm2,  T  =  500°K  and  V^c  3  lbl.6  V.  Table  9 
presents  the  results  of  the  large-signal  simulations  for  f  3  10  GHz 
and  T  3  500°K;  it  is  seen  that  by  increasing  E^q,  the  maximum  large- 
signal  efficiency  is  increased  to  over  29  percent.  The  v’  /V,  ratio 
is  now  0.8l  at  maximum  efficiency  compared  to  0.72  for  the  previous 
case  with  smaller  E  .  The  dc  operating  voltage  drops  approximately 

vO 

20  percent  from  VD„  3  0  to  V  3  92  V,  as  seen  in  Fig.  38b;  Fig.  38a 

nr  nr 

shows  that  the  hig’,. -efficiency  mode  has  become  more  critical  for  this 
structure,  that  is  the  efficiency  curve  is  more  peaked  and  the 
negative  conductance  curve  is  starting  to  increase  near  the  maximum 
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FIG.  37  (a)  HYBRID  DOPING  PROFILE  WHICH  FURTHER  INCREASES  E 
(b)  DC  SOLUTION  AT  T  =  500°K  AND  J,  =  1500  A/ cm2. 


Table  9 


Large-Signal  Results  for  the  Diode  of  Fig.  37a 


at 

T  =  500°K  and  f 

=  10 

GHz 

V 

RF  G 

B 

n 

Vdc 

Jdc 

p 

RF 

(V)  (U/cm1) 

(u/cm’) 

(%i 

lv) 

(A/cm1) 

(W/cm*) 

io  -56.712 

62.006 

1.34 

l4l 

1500 

2.8351x10* 

20  -44.932 

83.553 

4.27 

l4o 

1500 

8.9854x10* 

30  -35.161 

94.36 

7.61 

138 

1500 

15.821x10’ 

Uo  -28.585 

99.703 

11.05 

136 

1520 

22.867x10’ 

50  -23.558 

103.25 

14.56 

133 

1520 

29.446x10’ 

60  -19.52 

105.95 

18.06 

130 

1500 

35.135x10’ 

70  -16.522 

107.94 

21.43 

126 

1500 

40.478x10’ 

80  -lU.15 

109.41 

24.69 

121 

1520 

45.276x10’ 

85  -13.078 

110.2 

25.55 

118 

1570 

47.241x10’ 

87  -12.272 

110.86 

25.77 

117 

1550 

46.44x10’ 

88  -II.616 

111.47 

25.85 

116 

1500 

44.974x10’ 

89  -12.035 

110.98 

26.19 

115 

1580 

47.663x10’ 

90  -11.535 

111.5 

26.52 

114 

1540 

46.715x10’ 

91  -11. 1*9 

111.8.1 

29.00 

113 

1460 

47.572x10’ 

91  -11.69 

111.45 

28.10 

113 

1530 

48.401x10’ 

92  -11.832 

111.46 

29.08 

113 

1530 

50.071x10’ 

93  -10.689 

112.55 

28.11 

112 

1470 

46.224x10’ 

94  -  9.95U 

113.28 

26.52 

111 

1490 

43.974x10’ 

95  -  9.65 

113.56 

26,25 

111 

1490 

43.545x10’ 

96  -  8.7767 

114.17 

23.83 

110 

1540 

40.441x10’ 

98  -  5.8704 

116.65 

17.44 

109 

1490 

28.188x10’ 

99  -  1.4408 

120.12 

4.68 

105 

1430 

7.06x10’ 
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FIG.  38  LARGE-SIGNAL  RESULTS  FOR  THE  HYBRID  STRUCTURE  OF  FIG.  37a  AT  T  =  500°K,  f  =  10  GHz 
AND  J  =  1500  A/cm2. 


efficiency  point.  Therefore  it  appears  that  further  increase  of 
will  not  be  advantageous.  Figure  39  shows  the  diode  waveforms 
for  the  large-signal  solution  for  which  29.0-percent  efficiency  was 
obtained.  The  dip  in  the  current  near  270  degrees  is  brought  about 
by  the  inward  movement  of  the  depletion  edge  as  it  comes  to  meet  the 
electron  pulse.  -.Che  large  RF  voltage  modulation  possible  for  this 
structure  results  in  higher  efficiency  than  for  the  previous  hybrid 
double-drift  structures. 

The  effect  of  varying  the  uniform  doping  on  the  p-side  in 
hybrid  structures  is  examined  next.  Previously,  the  doping  wa3  set 
so  that  at  dc  the  electric  field  at  the  left-hand  boundary  was 
approximately  one- fourth  the  maximum  electric  field  value.  It 
is  not  clear  what  would  happen  if  the  doping  were  increased  resulting 
in  a  lower  field  at  the  left-hand  side.  The  avalanche  width  would 
be  decreased  which  should  be  beneficial;  on  the  other  hand  a  portion 
of  the  p-layer  may  become  undepleted  during  the  RF  cycle  which  would 
introduce  loss  due  to  the  low  mobility  for  holes. 

Figure  l*0a  shows  a  doping  profile  that  was  simulated  to  test 
the  effect  of  11^;  it  is  identical  with  the  profile  of  Fig.  3**a 
except  for  which  was  set  so  that  the  diode  is  just  punched  through 
at  dc,  as  shown  in  Fig.  40b.  For  the  structure  of  Fig.  3^a,  the 
maximum  efficiency  obtained  was  27.67  percent.  Figure  4l  shows  the 
waveforms  at  maximum  efficiency  for  the  just  punched  through  structure 
of  Fig.  **0a,  where  n  =  2k. 93  percent.  The  results  of  this  large- 
signal  simulation  show  that  much  of  the  p-side  is  undepleted  near 
270  degrees  in  the  cycle.  Therefore  it  is  desirable  to  design  the 
p-side  to  be  well  punched  through  at  dc. 
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FIG.  39  LARGE-SIGNAL  SOLUTION  FOR  THE  HYBRID  DIODE  OF  FIG.  37a 
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Simulations  were  also  carried  out  for  a  structure  with  lower 

than  that  of  Fig.  3**a  so  that  the  field  on  the  left-hand  side  was 

one-half  the  maximum  electric  field  at  dc.  The  maximum  large-signal 

efficiency  obtained  was  27.28  percent,  only  slightly  lower  than  for 

the  diode  of  Fig.  3ta;  however,  the  range  of  VRp  for  which  high 

efficiency  is  obtained  was  reduced.  These  results  indicate  that 

designing  so  that  the  field  on  the  left-hand  side  is  0.25  E  at 

max 

dc  is  desirable,  but  this  field  value  is  not  a  critical  factor. 

h. 3  Double-Read  Double-Drift  Structures.  The  double-d^ift 

diodes  considered  up  to  now  were  uniformly  doped  on  the  p-side  which 

resulted  in  a  large  avalanche  region  width.  Thi3  section  considers 

structures  with  high-low,  Read  doping  profiles  on  both  sides  of  the 

Junction.  Such  a  structure  can  theoretically  achieve  higher 

efficiencies  than  hybrid  double-drift  IMPATTs. 

Figure  lt2a  shown  the  first  double- Read  structure  that  was 

studied.  For  this  case  the  highly  doped  regions  on  both  sides  of 

the  junction  are  equal.  The  dc  solution  with  J^c  »  1250  A/cm2  and 

T  *  500°K  is  shown  in  Fig.  U2b .  The  structure  was  designed  so  that 

E.  (field  at  the  beginning  of  the  drift  region)  is  approximately 
to 

the  same  on  both  sides.  Also  the  electric  field  at  the  left-hand 
boundary  is  designed  above  the  value  at  the  right-hand  boundary  so 
that  the  diode  remains  punched  through  on  the  p-side  throughout  the 
entire  RF  cycle. 

Table  10  presents  results  of  large-signal  simulations  with 
f  *  10  GHz  and  T  =  500°K;  it  is  seen  that  efficiencies  of  over  36 
percent  are  obtainable.  Figure  h3  shows  plots  of  the  results  vs. 

RF  voltage  for  the  cases  with  Jdc  near  1250  A/ cm2.  Figure  U3a  shows 
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DOUBLE-READ  DOPING  PROFILE  AND  (b)  DC  SOLUTION  AT  T  =  500 °K  AID  =  1250  A/e*2 


Table  10 


Large-Signal  Results  for  the  Double-Read  Diode  of  Fig.  42a 
with  T  *  500°K  and  f  «  10  GHz 


VRF 

ill 

G 

(u/cm1) 

B 

(u/cm*) 

n 

i*> 

Vdc 

M 

Jdc 

(A/cm*) 

p 

RF 

(W/cm* ) 

10 

-58.003 

51.029 

2.09 

112 

1240 

2.9x10* 

20 

-U8.057 

88.187 

6.96 

111 

1240 

9.6112x10’ 

30 

-37.325 

101.32 

12.05 

no 

1260 

16.796x10* 

40 

-29.556 

107.9 

17.27 

109 

1260 

23.644x10’ 

50 

-23.779 

111.77 

22.55 

107 

1240 

29.723x10* 

6o 

-19.9 

113.85 

27.70 

104 

1240 

35.818x10’ 

65 

-I8.l4l 

111*.  71 

30. ll* 

102 

1250 

38.321x10’ 

70 

-16.222 

115.78 

32.37 

99.7 

1230 

39.741x10’ 

75 

-15.835 

115.68 

36.26 

96.7 

1270 

44.534x10’ 

76 

-14.782 

116.52' 

36.17 

95.8 

1230 

42.69x10’ 

77 

-14.89 

116.28 

35.6 

95.6 

1300 

44.138x10’ 

78 

-12.1*91 

118.58 

3l*.0l* 

94.3 

1180 

37.997x10’ 

80 

-10.382 

120.01 

29.  ll* 

92.0 

1240 

33.22x10’ 

82 

-65.84 

1240 

75 

-12.15 

119.05 

36.13 

95.1 

994 

34.17x10’ 

75 

-  8.6966 

122.13 

34.72 

93.6 

753 

24.458x10’ 

75 

-  U.9939 

125.3 

31.17 

92.1 

489 

14.045x10’ 

75 

-17.791* 

113.7 

33.57 

98.2 

1520 

50.044x10* 

76 

-19.93 

110.51* 

31.25 

98.8 

i860 

57.554x10* 

75 

-19.727 

111.1*6 

32.21 

99.2 

1740 

55-479x10’ 

75 

-18.539 

112.92 

33.08 

98.6 

1600 

52.138x10’ 

75 

-16.1*1*8 

115.15 

34.55 

97.6 

1370 

46.257x10’ 
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that  the  efficiency  and  negative  conductance  ourven  for  this  device 
are  very  well  behaved.  Figure  l)3b  shows  that  the  dc  operating  voltage 
decreases  approximately  l>t  percent  from  V  «  0  to  the  maximum 
efficiency  point.  Figure  UU  shown  the  computed  largo-signal  efficiency 
as  a  function  of  dc  current  density  for  conatunt  VD1,.  Tt  in  seen 
that  excellent  efficiency  ia  obtained  over  a  broad  range  of  current 
density.  Figure  h1)  shows  the  device  waveforms  for  the  most  efficient 
large-signal  simulation  with  n  =  3 (>.26  percent.  It  is  seen  that  the 
electron  pulse  is  efficiently  collected  resulting  in  a  very  small 
current  at  90  degrees  in  the  cycle.  Figure  H6  shows  the  electric 
field  profile  and  particle  concentrations  at  two  times  in  the  cycle; 
Fig.  Ii6b  is  Just  before  the  pulse  i:;  collected.  The  p-side  is 
always  punched  through  during  the  cycle  due  to  the  space  charge  of 
the  generated  hole  pulse.  The  generated  electron  pulse  is  larger  than 
the  background  doping,  as  it  should  bo  for  efficient  precollection. 
Figure  l*6a  shows  that  the  device  is  nearly  drained  out  by  the  end 
of  the  cycle. 

In  the  next  structure  of  Fig.  hja.,  E  was  reduced  to 
1.8  x  105  V/cm  on  both  sides  so  that  the  diode  is  not  punched  through 
on  the  n-side,  as  shown  by  the  dc  solution  of  Fig.  UTb.  Based  on  the 
hybrid  simulation  results,  it  is  expected  that  decreasing  will 
decrease  the  efficiency  of  the  device.  The  simulation  results  shown 
in  Table  11  confirm  that  the  efficiency  drops  from  over  36  percent 
for  the  punch- through  case  to  26.^1  percent  for  the  nonpunch-through 
case.  Figure  U8  shows  the  waveforms  for  the  V  =  U9  V  solution; 

r\r 

the  current  at  90  degrees  is  considerably  larger  than  for  Fig.  1+5 
because  the  n-side  is  longer  and  the  pulse  is  not  entirely  collected 
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DOUBLE  DRIFT 
VRF=75. 

JDC-1270.  A/CM*k2 
FREQUENCY  =  10.0  GHZ 


PHASE  ANGLE,  DEGREES 


PHASE  ANGLE,  DEGREES 

FIG.  i»5  LARGE-SIGNAL  SOLUTION  FOR  THE  DOUBLE- READ  STRUCTURE 
OF  FIG.  l»2a  WITH  T  =  500°K  AND  f  =  10  GHz  RESULTING 
IN  36.26-PERCENT  EFFICIENCY. 
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Table  11 


Large-Signal  Results  for  the  Double-Read  Structure  of  Fig.  47a 
at  T  ■  500°K  and  f  »  10  GHz  • 


v 

RF  G 

B 

n 

Vdc 

^dc 

PBF 

(V)  {0/  cm*) 

(u/cm1) 

ill 

M 

(A/cm‘) 

Wcm») 

5.0  -51.758 

45.874 

0.68 

76.5 

1250 

646.93 

uo.o  -21,252 

106.38 

18.67 

73.9 

1230 

17.001x10* 

1*2.0  -20.527 

106.7 

19.68 

73.8 

1250 

18.105x10* 

41*.  0  -19.544 

107.64 

20.58 

73.6 

1250 

18.918x10* 

46.0  -19.233 

108.85 

22.46 

77.0 

1240 

20.348x10’ 

48.0  -20.151 

111.04 

26.02 

71.7 

1240 

23.214x10* 

49.0  -19.637 

111.48 

26.41 

71.6 

1250 

23.575x10* 

50.0  -18.808 

112.12 

26.45 

71.6 

1240 

23.51x10* 

DOUBLE  DRIFT 
VRF=49. 

JDO1250.  R/CMmm2 
FREQUENCY  =  10.0  GHZ 


PHASE  ANGLE,  DEGREES 


PHASE  ANGLE,  DEGREES 


FIG.  U8  LARGE-SIGNAL  SOLUTION  FOR  THE  DOUBLE-READ  STRUCTURE 

OF  FIG.  UTa  WITH  T  =  500°K  RESULTING  IN  26 . Ul-PERCENT 
EFFICIENCY. 
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at  the  end  of  the  previous  cycle.  Also  the  VRp/Vdc  ratio  is  smaller 
than  for  the  punch-through  case  which  contributes  to  lower  efficiency. 
The  n-region  can  be  shortened  by  decreasing  E  further  which  would 

wO 

reduce  the  current  during  the  first  half  of  the  cycle;  however,  the 
VRp/Vdc  ratio  would  also  be  further  reduced  and  the  efficiency  would 
still  be  well  below  that  obtainable  from  the  punched- through  device. 

2*  Incorporation  of  Dynamic  Temperature  Effects  in  IMP ATT  Simulations 
Up  to  now  it  has  been  assumed  that  the  diode  temperature 
remains  constant  as  VRp,  Jdc  and  f  are  varied  in  the  large-signal 
IMPATT  simulations.  Thin  is  not  quite  true  since  the  diode  temperature 
is  given  by 

T  =  300  ♦  <1  -  n>VdcJdcKth  ,  (00) 

where  T  is  the  temperature  in  °K,  n  is  the  fractional  efficiency, 

V^c  is  the  dc  operating  voltage,  is  the  dc  operating  current  in 
A/cm2,  and  R  ^  is  the  diode  thermal  resistance  in  (°K-cm2)/W. 

It  is  seen  that  as  the  efficiency  and  dc  parameters  of  the  large- 
signal  solution  vary,  the  diode  temperature  will  also  change. 

That  this  effect  should  be  included  in  the  IMPATT  simulation 
became  apparent  when  a  series  of  large-signal  runs  were  carried  out 
at  constant  dc  operating  voltage.  The  structure  of  Fig.  21a  was 
previously  simulated  for  varying  between  900  and  1000  A/cm2,  as 
shown  in  Table  6.  At  the  maximum  efficiency  point,  the  dc  operating 
voltage  was  61.9  V.  The  large-signal  IMPATT  program  was  modified  to 
run  at  constant  dc  voltage  and  this  structure  was  again  simulated  for 
different  RF  voltages  but  with  V^c  =  constant  =  6l.9  V.  Therefore  in 
this  mode  of  operation,  J^c  adjusts  itself  to  whatever  value  is 


consistent  with  the  given  voltage  drive.  Table  12  shows  the  results 
obtained  and  Fig.  1*9  shows  a  plot  of  efficiency  and  negative 
conductance  vs.  V  It  is  seen  that  the  diode  does  not  begin  to 
conduct  until  Vpp  =  40  V  is  reached.  However  all  the  simulations 
in  Table  12  were  run  at  a  constant  temperature  of  T  *  500°K;  this 
is  obviously  not  consistent  with  Eq.  20  since  below  VRF  -  kO  V  the 
device  is  not  dissipating  power  and  should  therefore  be  at  room 
temperature . 

At  this  point,  the  program  was  modified  as  follows.  The 
capability  of  reloading  the  diode  material  parameters  was  incorporated. 
The  A  and  b  coefficients  in  the  ionization  rate  expression  given  in 
Table  1  were  assumed  to  vary  linearly  with  temperature  between  the 
values  given  at  300°K  and  500°K  based  on  the  results  of  Hall  and 
Leek.8  The  low- field  mobilities  were  assumed  to  vary  as 

n(T)  «  uo(300/T)n  ,  (21) 

where  uq  is  the  value  at  300°K  in  Table  1,  T  is  the  temperature  in  °K 
and  n  is  an  exponent  chosen  so  that  p(500°K)  agrees  with  the  value  in 
Table  1.  All  the  other  material  parameters  in  Table  1  were  assumed 
to  vary  linearly  with  T  such  that  the  values  at  T  *  300°K  and  T  *  500°K 
agree  with  the  values  given  in  Table  1. 

The  algorithm  may  be  described  a3  follows.  At  the  start  of 
the  simulation  a  material  parameter  temperature  T^  and  thermal 
resistance  R^  are  assumed  and  an  RF  solution  is  generated  at 
yielding  values  for  n»  V^c  and  J^.  T  is  then  calculated  according 
to  Eq.  20;  in  general,  T  j  T^.  At  this  point  the  material  parameters 
are  reloaded  at  a  new  temperature  T^  and  a  new  RF  solution  at  T^ 
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Table  12 


Large-Signal  Results  for  the  Diode  of  Fig.  21a  Using  a  Constant 
Dc  Voltage  of  6l.9  V  and  T  ■  500°K 


v 

RF  G 

(V)  ( U/cm1) 

B 

(U/CK1) 

n  (%) 

f 

(GHz) 

Jdc 

(A/cm1) 

P op 

(W/c  «•). 

1*0.0  ~o.o 

10.0 

-0.0 

1*5.0  -  2.255 

217.1*9 

16.35 

10.0 

226 

2.2832x103 

1*7.0  -  l*.l6ll* 

217.09 

21.1*7 

10.0 

31*6 

1*.  5962x103 

1*9.0  -  7.1*25 

215.1* 

27.13 

10.0 

531 

8.9136x103 

51.0  — ll* .  86l 

210.0 

3U.80 

10.0 

897 

19.326x103 

51.0  -17.623 

196.01* 

38'.  3l* 

9.5 

966 

22.919x103 

52.0  9.1273 

171.55 

-  1* .  »*3 

10.0 

1*500 

-12.3l*xl03 

53.0  30.809 

135.1*5 

-11.01 

10.0 

6350 

-1*3.27x103 

FIG.  49  LARGE-SIGNAL  RESULTS  FOR  THE  DIODE  OF  FIG.  2ia 
AT  T  «  500°K,  f  »  10  GHz  AND  V  =  6l.9  V. 
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PERC 


yields  n’»  V!^  and  J!^.  T'  is  calculated  next  as  follows; 


*’  ■  300  +  (i  -  n '  ^Vd0J4eBtlj  • 


(22) 


where  in  general  T*  4  T^.  Now  from  the  parameter 


AT  T*  -  T 

AT  T'  -  T 

P  P  P 


(23) 


a  measure  is  obtained  of  how  the  temperature  according  to  Eq.  20 
varies  as  a  function  of  device  temperature  T^.  If  T  i3  assumed  to 
vary  linearly  with  T^,  the  following  expression  is  obtained: 


T  (T  ) 
v  pv 


T  -  T' 


T  -  T* 

P  P 


rptrp  _  ipipf 

...  + - E _ E 

pv  T  -  T' 

P  P 


(2k) 


where  Ty  and  T^y  are  variables  representing  the  temperature 
calculated  according  to  Eq.  20  and  the  material  parameter  temperature, 
respectively.  If  Eq.  is  set  equal  to  T^v>  a  new  material 

parameter  temperature  T^  is  obtained  which  is  a  guess  (assuming 
linearity  in  the  T  vs.  T  relation)  of  the  required  device  temperature: 


fpH 

P 


TT» 

_ R 

T  -  T’ 


-  T'T 


_ E _ 

-  T  +  T*  * 

P  P 


(25) 


The  material  parameters  are  reloaded  at  T^  and  a  third  BP  solution  is 
calculated  yielding  n"»  Vjjc  and  J'^.  The  following  is  then  calculated: 

T"  =  300  +  (1  -  n")VdCJ£cRth  •  (26) 

Most  of  the  time  it  is  found  that  T"  ~  T"  to  within  an  acceptable 

P 

tolerance  and  the  simulation  is  terminated.  If  T"  4  T”,  the  process  is 
continued  until  a  thermally  consistent  solution  is  obtained. 
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For  example,  the  structure  of  Fig,  15a  is  considered  which  was 

previously  run  at  constant  T  *  500°K  and  constant  J,  s  1200  A/cm2. 

dc 

The  best  efficiency  (see  Table  5)  was  obtained  for  -  66.1  V 

and  V__  ■  51,9  V.  A  series  of  simulations  with  constant  V.  »  66.1  V 
nr  aC 

was  carried  out  for  this  structure  with  variable  temperature.  The 
diode  thermal  resistance  was  chosen  to  be  consistent  with  the  optimal 
solution  just  mentioned,  i.e.. 


so  that  the  solution  for  VRR  »  51.9  V  should  yield  T  »  500°K  and 
n  ■  3^.18  percent  as  before. 

Table  13  shows  the  results  up  to  the  maximum  efficiency  point*, 

Fig.  50  shows  plots  of  these  results  at  f  ■  10  GHz  which  show 

that  the  device  heats  up  from  approximately  1*50°K  at  dc  to  500°K  at 

maximum  efficiency.  These  curves  are  In  sharp  contrast  with  those  of 

Fig.  1*9  where  the  temperature  was  constant  at  T  »  500°K.  It  is  seen 

that  entirely  erroneous  results  may  be  obtained  if  the  temperature 

variation  is  not  included.  From  Fig.  50a  it  is  seen  that  operating 

the  diode  at  constant  considerably  flatt*-r.3  out  the  negative 

conductance  curve  which  would  be  advantageous  for  amplifier  applications. 

However,  as  shown  in  Fig.  50b  the  dc  current  increases  very  rapidly 

near  the  maximum  efficiency  point  so  that  probably  some  sort  of 

current  limiter  would  be  required  to  prevent  burnout. 

The  structure  of  Fig.  15a  was  next  simulated  at  constant 

J.  *  1200  A/cm2  but  with  variable  temperature.  The  diode  thermal 
dc 

resistance  cannot  be  chosen  consistent  with  the  maximum  efficiency 


--n y ' -.'Vr.':'.y.'^ 5,V*S'-T^?,v  ,  *;c 


Table  13 

Large-Signal  Results  for  the  Diode  of  Fig.  15a 

with  f  «  10  0Hz,  V(lc  *  66.1  V  and 

R..  ■  3.963  x  10“  3  (°K-cm2/W) 
tn 


v 

RF  G 

M  (o/cm*) 

B 

(tl/cm1)  n  (%) 

Jdc 

(A/cm* )  T(°K) 

(W/cm*) 

5.0  -56.895 

183.1 

1.76 

610 

457.7 

0.7111x10* 

10.0  -44.336 

189.92 

5.2l* 

640 

456.4 

2.2167x10* 

15.0  -33.617 

191).  95 

8.93 

641 

456.4 

3.7819x10* 

20.0  -28.13 

197.0 

12.63 

674 

455.0 

5.626x10' 

25.0  -23.895 

198.29 

16.3 

693 

455.0 

7.467x10’ 

30.0  -21.382 

198.9 

19.65 

741 

456.5 

9.622x10’ 

1*0.0  -17.59 

200.1*3 

2U.97 

852 

468.7 

14.072x10’ 

1*5.0  -16.938 

200.96 

27.62 

939 

477.5 

17.149x10’ 

1*8.0  -16.255 

201.69 

28.53 

993 

484.7 

18.725x10* 

50.0  -16.017 

202.05 

29.38 

1030 

491.9 

20.021x10’ 

51.9  -19.1)2 

199.92 

3l*.l8 

1160 

500.0 

26.155x10* 

-98- 


solution  as  before,  since  for  constant  current  bias  the  diode  heats 
up  as  V__  is  reduced  so  that  it  would  bum  out  using  the  R..  for 

nr  UR 

the  simulations  of  Table  13.  (The  diode  heats  up  because  the 
efficiency  in  Eq.  20  decreases  and  there  is  still  significant  dc 
power  flowing.)  Instead,  the  thermal  resistance  was  chosen  consistent 
with  the  dc  solution  at  T  *  500°K  so  that 


"th 


500  -  300 
( 1200 )( 77.36) 


2.1544  x  10" 3  (°K-cm2/W) 


(28) 


Table  14  shows  the  simulation  results  and  Fig,  51  presents 

plots  of  these  results.  These  curves  should  be  compared  with  those 

of  Fig.  17  for  the  3a»e  J.  but  constant  T  »  500°K.  The  maximum 

dc 

efficiency  is  approximately  the  same,  but  the  efficiency  curve  of 
Fie.  51a  is  much  better  behaved  than  that  of  Fig.  17a  due  to  the 
temperature  variation.  Ficure  ,51b  shows  that  the  device  cools  down 
with  increasing  VDT1  when  constant  current  bias  is  used  since  the 
efficiency  is  increasing.  For  high-efficiency  GaAs  IMPATTs,  it  is 
quite  important  to  design  the  thermal  resistance  properly.  Since 
V^c  significantly  increases  with  decreasing  VRp  and  n  significantly 
decreases,  Eq.  20  predicts  that  the  diode  can  easily  bum  out  below 
the  optimum  value  if  R^h  is  too  large.  For  example,  this  diode 
bums  out  using  constant  current  bias  if  the  R^  used  for  the  constant 
voltage  bias  case,  3.96  x  10"3  (°K-cm2/W),  is  used  instead  of 
2.154  x  10’3  (°K-cm2/W). 


i 
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Table  l4 


Large-Signal.  Results  for  the  Structure  of  Fig.  15a 


with  f  = 

:  10  GHz 

*  ^dc 

1200  A/cm2 

c 

anclRth 

*  2.154 

x  10“  3 

(°ft-cra2/W) 

v 

RF  G 

(V)  (U/cm*) 

B 

(u/cm*) 

n  (%) 

Jd^ 

(A/cm*) 

Vdc 

(V) 

/ 

T  (°K)  ' 

P 

RF 

(W/cm*) 

5.0  -103.67 

155.31 

1.1*1 

1190 

77.2 

500 

1. 296x10* 

10.0  -  73.839 

179.59 

4.11 

1190 

75-5 

490 

3.692x10’ 

15.0  -  59.051* 

184.85 

7.72 

1180 

72.7 

473.3 

6.643x10* 

20.0  -  1*9.955 

187.52 

11.51 

1220 

71.0 

461.9 

9.991x10’ 

25.0  -  1*2.162 

190.08 

15.82' 

1200 

69.1 

451.9 

13.175x10’ 

30.0  -  37.902 

189.1*1* 

21.23 

1210 

66.1 

435.33 

17.056x10* 

35.0  -  34.311 

187.39 

27.1*5 

1240 

61.8 

417.35 

21.015x10’ 

36.0  -  32.99 

187.25 

28. 94 

1220 

60.4 

412.35 

21.377x10*. 

37.0  -  31.389 

187.1*5 

30.85 

1180 

58.9 

407.35 

21.486x10’ 

38.0  -  32.225 

185.19 

33.55 

1220 

57.0 

400.35 

23.266x10*, 

39.0  -  29.919 

187.68 

3»*.09 

1180 

56.6 

400.35 

22.753x10’ 

1*0.0  -  27.738 

189.83 

33.49 

1200 

55.2 

395 

22.19x10’ 

1*1.0  -  21 .21 

192.99 

32.64 

1190 

54.6 

395 

21.188x10’ 

1*2.0  -  22.587 

196.51 

30.91 

1190 

54.2 

395 

19.932x10’ 

1*3.0  -  21.173 

198.1*8 

30.17 

1190 

54.7 

400 

19.574x10’ 
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6.  Effects  of  Material  Parameters 

The  simulations  presented  in  this  report  so  far  have  used 
•the  ionization  rates  determined  by  Hall  and  Leek8  and  the  particle 
velocities  from  the  references  presented  in  Bauhahn's  report.6 
Recently,  new  measurements  of  ionization  rates  and  hole  velocity  in 
GaAs  were  published  by  Raytheon.9  Table  15  shows  the  new  material 
parameters  based  on  the  Raytheon  results,  If  the  ionization  rate 
expression  is  evaluated  and  compared  with  the  values  given  in  Table  1, 
it  is  found  that  the  ionization  rates  published  by  Raytheon  are 

«  t  ' 

slightly  lower  than  the  Hall  and  Leek  rates.  For  example,  at  T  =  300°K, 
E  *  1000  kV/cm,  the  Hall  and  Leek  ionization  rate  is  1.1*78  x  10 5  cm”1, 
whereas  the  Raytheon  ionization  rate  is  1.202  x  105  cm”1.  To  study 
the  effect  of  the  lower  ionization  rat03,  the  values  in  Table  15  were 
incorporated  in  the  simulation  and  a  one-sided,  high-low  diode  very 
similar  to  the  structure  of  Fig.  6  was  simulated.  Previously,  a 
maximum  efficiency  of  27.9  percent  was  obtained  using  the  Hall  and  Leek 
rates;  using  the  rates  reported  by  Raytheon,  a  maximum  efficiency  of 
27.25  percent  was  obtained.  This  comparison  indicated  that  the 
reduction  in  ionization  rates  has  a  negligible  effect  on  the  calculated 
results . 

More  significant  is  the  lower  hole  saturated  velocity  in  Table  15. 
Of  course,  for  one-sided  diodes,  incorporation  of  these  values  has 
negligible  effect  on  the  results.  However,  for  double-drift  diodes  the 
length  of  the  p-side  must  be  reduced  in  order  to  obtain  similar  results. 
Comparison  with  Table  1  shows  that  the  widths  of  all  the  p-layers 
simulated  so  far  would  be  reduced  by  a  factor  of  approximately  0.8 
if  Raytheon's  saturated  hole  velocity  is  incorporated. 


7.  Power  Output  from  IMPATT  Devices 

The  simulations  in  this  report  were  carried  out  for  normalized 
(per  unit  area)  cases,  so  that  absolute  impedance  ar.d  power  levels 
have  not  been  specified.  In  this  section,  achievable  thermal  resistances 
for  single-  and  double-drift  IMPATTs  are  used  to  estimate  the  power 
output  from  various  devices  already  studied.  It  is  important  to  note 
that  the  results  presented  here  are  at  the  maximum  efficiency  points 
for  the  various  structures.  Higher  powe~  outputs  (i.e.,  under 
pulsed  conditions)  can  be  achieved  if  very  short  pulses  are  utilized 
and  where  device  heating  is  not  a  problem. 

Table  16  presents  a  summary  of  the  cases  considered  in  this 
section.  In  the  diode  type  column,  SD  denotes  single  drift  and  DD 
denotes  double  drift;  PT  and  NPT  denote  punch  through  and  nonpunch 
through,  respectively.  DR  denotes  double  Read  and  HYB  denotes  hybrid. 

For  each  structure,  large-sighal  results  for  the  maximum  efficiency 
point  are  given.  In  Table  .‘.6,  the  diode  area  is  chosen  such  that  the 
negative  resistance  exhibited  by  the  diode  is  -  1  ft  for  each  case;  this 
is  because  we  are  assuming  a  circuit  load  resistance  of  1  ft  which  is 
reasonable  in  practice.  Hence  the  output  powers  given  in  Table  l6 
are  the  maximum  which  could  be  achieved  for  this  value  of  load 
resistance  if  the  required  thermal  resistance  can  be  obtained.  6 

n 

is  the  required  thermal  resistance  for  CW  operation  so  that  the 
device  temperature  does  not  exceed  200°C. 

For  the  one-sided,  uniformly  doped  IMPATT  structure  of  Fig.  4, 
an  area  of  2.73  x  10”"  cm*  corresponds  to  a  device  negative  resistance 
of  -  1  ft .  If  the  external  circuit  can  match  the  diode  at  this 
impedance  level,  the  resulting  CW  RF  power  is  2.28  W  and  the  diode 


thermal  resistance  required  for  operation  at  500°K  is  15.95  (°K/W) . 
Examination  of  the  thermal  resistance  curve  in  Szel#  shows  that  this 
value  of  thermal  resistance  can  be  achieved. 

For  the  one-sided,  high-low,  nonpunch- through  diode  for  which 
large-signal  results  are  presented  in  Table  2,  the  most  efficient 
solution  was  with  VRp  a  29.0  V,  for  which  n  =  27.9  percent.  Again, 
matching  to  1-fl  resistance,  the  resulting  values  are  Pe_  =  13.5  W  (CW), 

ni* 

A  «  8.59  x  10”"  cm1,  Rth  =  5«73°K/W.  According  to  Sze,  this  value  of 
thermal  resistance  is  reasonable. 

Next,  the  one-sided,  low-high-low,  punch-through  diode  for 
which  large-signal  results  are  presented  in  Table  6  is  considered. 

Going  through  the  same  procedure  yields  PRfl  =  10.2  W  (CW)  and  » 
12.18°K/W  at  the  most  efficient  point;  this  value  of  RTH  is  reasonable 
according  to  Sze’s  curve.  It  is  noted  that  the  predicted  CW  power  is 
less  than  for  the  previous  high-low,  nonpunch-through  diode;  this  may 
counter  one's  expectation,  since  the  efficiency  for  the  punch-through 
diode  is  considerably  higher.  The  discrepancy  is  explained  by  noting 
that  the  current  density  for  the  simulation  of  the  nonpunch-through 
structure  was  J^c  =  1280  A/cm’,  while  for  simulation  of  the  punch- 
through  structure  it  was  =  957  A/cm*.  If  the  punch-through 
structure  were  simulated  at  a  higher  J^,  it  would  be  capable  of 
generating  as  much  power  as  the  nonpunch- through  structure;  also  the 
R,pjj  requirement  would  be  less  stringent  (higher)  for  the  punch-through 
structure  since  the  efficiency  is  higher  (dissipated  power  is  lower). 

The  simulations  of  Table  16  were  not  carried  out  with  the  goal 
in  mind  of  evaluating  the  maximum  power  generating  capabilities  of 


each  structure,  but  rather  of  predicting  maximum  efficiencies  obtainable 
at  10  GHz.  Therefore  a  direct  comparison  of  power  levels  is  not  always 
meaningful,  as  seen  in  the  previous  paragraph.  The  results  in  Table  16 
do  indicate  the  power  expected  for  each  structure  under  the  particular 
conditions  of  the  simulation. 

Next,  the  power  output  from  double-drift  IMPATTs  is  considered. 
These  diodes  cannot  typically  be  matched  to  1-n  resistance  because  the 
required  thermal  resistance  is  unrealistically  low  for  CW  operation. 

The  thermal  resistance  for  double-drift  diodes  is  given  approximately 


In  T  4t  T 

rth  =  STd  *  Rpkg  +  — —  +  — ^  »  (29) 

TH  *Khad  PKG  120ffdi  300ird* 

where  d  is  the  device  diameter  in  cm,  Khg  is  the  thermal  conductivity 
of  the  heat  sink  material  in  W/cm-°K,  Rp^,  is  the  portion  of  the  thermal 
resistance  due  to  the  package  and  bonding  (°K/W),  l  is  the  buffer  or 

i 

substrate  thjLckness  (cm),  T  is  the  average  temperature  (°K)  of  the 

i 

buffer  layer,  l  is  the  active  layer  thickness  (cm)  between  the 

i 

junction  and  the  buffer  layer,  and  T  is  the  average  temperature  (°K) 


of  the  active  layer. 


For  these  calculations,  the  thermal  conductivity  of  a  diamond 

heat  sink  is  taken  as  11.7  W/cm-°K  and  that  of  a  copper  heat  sink 

is  3*9  W/cm-°K.“  Also,  Rn,,„  is  assumed  to  have  the  constant  value  of 

rl\u 

2.1  °K/W.“  It  is  assumed  throughout  that  the  buffer  layer  thickness 
l  -  10-'  cm  and  that  the  average  buffer  layer  tenperature  is  T  =  450°K. 

t  » 

The  active  layer  is  assumed  to  be  the  p-layer  and  the  thickness  l  is 

2 

reduced  by  the  factor  0.8  from  the  values  specified  in  the  figures, 
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consistent  with  the  recent  measurement’ of  saturated  hole  velocity  in 
GaAs  discussed  in  Section  6.  The  average  temperature  in  the  active 
layer  is  assumed  to  be  T  =  500°K. 

i 

The  first  term  in  Eq.  29  is  the  spreading  term,  which  occurs 
at  the  interface  between  the  diode  and  heat  sink.  This  component  of 
the  thermal  resistance  can  be  reduced  by  using  multiple  mesa  arrays 
or  ring  geometries.  “  Figure  52  illustrates  the  reduction  of  the 
spreading  term  which  is  possible  by  using  a  ring  geometry.  “  In  these 
calculations  if  a  ring  structure  is  used,  it  is  assumed  to  reduce  the 
spreading  term  by  the  factor  0.55,  which  is  readily  achievable  accord¬ 
ing  to  Fig,  52. 

For  example,  if  the  active  region  length  t  =  2.5  x  10"”  cm 

i 

is  chosen,  Fig.  53  presents  R^  as  a  function  of  device  diameter  for 
different  heat  sink  materials  and  device  geometries.  First,  the 
uniformly  doped,  nonpunch-through  double-drift  structure  is  considered 
for  which  large-signal  results  are  given  in  Table  7.  The  most  efficient 
solu+ion  obtained  at  10  GHz  is  n  =  15*47  percent,  Vdc  =  117  V  and 
J^c  =  984  A/cm1.  The  reduced  p-layer  thickness  for  this  structure  is 
l  =  3*0  x  KT’  cm.  Since  the  operating  temperature  is  500°K, 

t 

500  =  300  +  (1  -  0.1547)(117)(984)it  (dV4)RTH  .  (30) 

It  is  noted  from  Table  16  that  if  the  device  is  matched  to  1  8  circuit 
resistance,  a  thermal  resistance  of  3*63°K/W  is  required  with  a  diode 
diameter  of  10.57  mils;  however,  Fig.  53  shows  that  at  this  diameter, 
a  thermal  resistance  of  3»63°K/W  is  not  achievable.  To  find  the 
expected  power  output,  Eqs.  29  and  30,  which  constitute  two  equations 
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FIG.  53  R^n  FOR  DOUBLE-DRIFT  G&As  DIODES  VS.  DIAMETER  FOR 
(1)  COPPER  HEAT  SINK,  CIRCULAR  GEOMETRY;  (2)  COPPER 
HEAT  SINK,  RING  GEOMETRY;  (3)  DIAMOND  HEAT  SINK, 
CIRCULAR  GEOMETRY;  AND  (4)  DIAMOND  HEAT  SINK,  RING 
GEOMETRY.  <1  =  2.5  x  KT"  cm  AND  l  =  lO-*  cm) 
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for  the  two  unknowns  R^  and  d,  are  solved  simultaneously.  The 
maximum  expected  power  is  obtained  with  the  equality  sign  in  Eq.  30. 

Table  IT  shows  the  maximum  expected  RF  power  for  the  double¬ 
drift  structures  for  different  heat  sink  materials  and  mesa  geometries. 
For  the  uniform  double-drift  diode  under  the  most  favorable  conditions, 
5.91  W  CW  output  power  is  expected,  or  nearly  twice  the  power  achievable 
with  the  uniform,  single-drift  diode. 

Next,  the  nonpunch-through,  hybrid  double-drift  diode  of  Fig.  28a 
is  considered.  The  solution  for  this  device  yielded  n  3  18.72  percent, 

=  104  V,  and  Jdc  =  1470  A/cm* .  The  reduced  p-layer  thickness 
used  in  Eq.  29  is  1  =  2.4  x  lO"*  cm.  The  equation  corresponding  to 

i 

Eq.  31  for  this  device  is: 


500  =  300  +  (1  -  0.l872)(104)(l470Md*/4)RTH  .  (31) 


Simultaneous  solution  of  Eqs.  29  and  31  yields  the  maximum  expected 
generated  powers  shown  in  Table  17. 

From  these  results,  it  appears  that  there  is  hardly  any 
advantage  to  using  a  hybrid  nonpunch-through  structure  instead  of  a 
uniform  structure.  However,  it  is  again  stressed  that  these  two 
structures  were  driven  with  considerably  different  dc  current  densities, 
and  the  goal  was  to  find  the  maximum  efficiency  point  rather  than  the 
maximum  power  point.  It  is  quite  possible  that  using  a  lower  J^c 
for  the  double-drift  hybrid  simulation  would  result  in  higher  expected 
PDT?;  this  is  explained  with  reference  to  Eq.  31*  If  Is  decreased 
by  some  factor  and  everything  else  remains  the  same,  the  device  area 
could  be  increased  by  the  inverse  factor,  which  would  result  in  the 
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same  PRp.  However  in  the  region  of  Fig.  53  where  R^  decreases 
rapidly  with  diameter,  the  area  could  be  increased  much  more  than  by 
the  inverse  factor,  resulting  in  a  greater  Ppp.  Of  course,  the  situation 
is  complicated  by  the  fact  that  n,  V^,  etc.,  change  as  Jdc  is 
varied.  The  important  point  is  that  the  values  given  in  Table  17 
are  not  to  be  interpreted  as  the  maximum  possible  CW  power  points, 
but  rather  as  typical  values  obtained  from  simulations  near  the 
maximum  efficiency  point. 

A  similar  procedure  was  carried  out  for  three  punch -through 
double-drift  hybrid  structures  at  different  levels  of  punch  through — the 
slightly  punched  through  hybrid  of  Fig.  33&>  for  which  E^q  = 

1.7  x  10s  V/cm;  the  punch-through  hybrid  of  Fig.  3^a>  for  which 
EfQ  =  2.23  x  10s  V/cm;  and  the  strongly  punch-through  hybrid  of 
Fig.  37a,  with  EtQ  =  2.68  x  10s  V/em.  From  the  results  given  in 
Table  17,  it  is  seen  that  it  is  advantageous  to  design  the  hybrid 
to  be  punched  through.  In  fact,  nearly  twice  the  power  of  the 
nonpunch-through  hybrid  can  be  obtained  by  the  punch-through  hybrid. 

Next,  two  double-Read  structures  were  examined— the  punch- 
through  structure  of  Fig.  l*2a  and  the  nonpunch-through  structure  of 
Fig.  U7a.  It  is  seen  that  the  punch-through  double-Read  structure 
is  capable  of  generating  more  power  than  any  other  structure  considered 
at  10  GHz.  It  may  be  surprising  to  compare  the  generated  power  for 
the  nonpunch-through  double  Read  with  that  of  the  punch-through 
hybrid.  The  efficiencies  are  about  equal,  but  the  double  Read  can 
generate  significantly  more  RF  power.  This  is  explained  by  noting 
from  Table  16  that  V^c  =  107  V  for  the  punch- through  hybrid  and 
Vdc  =  71.6  V  for  the  nonpunch-through  double-Read  structure.  Since 


-114- 


Vdc  is  significantly  smaller  for  the  double  Read,  it  can  operate 
with  a  larger  area  than  the  punch-through  hybrid.  Since 
rapidly  decreases  with  area,  the  ratio  of  the  device  areas  is  much 
larger  than  the  ratio;  hence,  the  double  Read  can  operate  at 
higher  power  levels. 

8.  Summary  and  Conclusions 

The  properties  of  various  types  of  doping  profiles  for  both 
single-  and  double-drift  GaAs  IMPATT  diodes  operating  near  10  GHz 
have  been  investigated.  From  the  results  given  in  Table  17,  the 
structures  may  be  rated  as  follows,  in  order  of  increasing  power 
generating  capability  and  efficiency:  single-drift  uniform  nonpunch- 
through,  single-drift  Read  nonpunch-through,  single-drift  Read 
punch-through,  double-drift  uniform  nonpunch-through,  double-drift 
hybrid  nonpunch-through,  double-drift  hybrid  punch- through,  double- 
Read  nonpunch- through,  and  double-Read  punch- through.  It  is  hoped 
that  these  results  will  provide  proper  guidelines  for  designing 
such  devices.  Tables  16  and  17  provide  a  summary  of  the  power  levels 
which  can  be  expected  from  the  various  structures  considered  here. 
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